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Fundamental Aspects of Solids‘(Gas Flow 


Part I: Introductory Concepts and Idealized Sphere Motion 
in Viscous Regime’ 


L. B. TOROBIN? and W. H. GAUVIN? 


This article introduces a series of analytical and 
critical literature surveys which deal first with the 
individual fundamental momentum transfer phenom- 
ena associated with particulate motion as affected by 
particle shape, roughness, acceleration, concentration 
and rotation, as well as the turbulence characteristics 
and boundary conditions of the transporting fluid. 
This information will then be used as a basis of 
discussion and interpretation of the many investiga- 
tions of multiparticle solids-gas conveying systems 
and contacting techniques reported in the literature. 

This first paper reviews briefly some pertinent 
basic concepts of fluid mechanics, and then considers 
the development of the fluid velocity field and drag 
for single spherical particles moving rectilinearly at 
pre-wake-forming Reynolds Numbers through an in- 
finite and undisturbed fluid. 


O' all the physical phenomena which, in various ways, 
affect man’s daily life, his environment and _ his 
activities, few exert as profound an influence as those 
associated with the occurrence of solids-gas systems. 
Whether evidenced in the form of large-scale meteoro- 
logical disturbances, such as sand and snow storms, or on 
a more modest scale, such as the dispersion of soot 
particles from a factory chimney, or finally in such short- 
range action as the entrainment of dust by vacuum- 
cleaning equipment, the behavior of solid particles in the 
presence of a gas stream significantly affects human 
welfare and economy. The range and variety of processes 
in which such systems play an important, if not a 
dominating part, is truly astounding: the spraying of 
insecticides, the dissemination of pollen, the motion of 
aircraft, rockets and missiles, the scattering of radioactive 
fallout, the inhalation of toxic dusts, the shifting of sand 
dunes, the erosion of soils, the color of the sky at tw ilight 
~ these are but a few random examples of their occurrence. 


The earliest harnessing of these systems for industrial 
uses was associated with materials handling and pneumatic 
Industrial equipment for this was 


convey ing. purpose 
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probably first reported in the literature in 1897 by 
Duckam ‘)), although it _ wee employed in small- scale 
operations as early as 1866 ‘?). Whereas the installations 
of the beginning of the century were used primarily for 
low-speed vacuum movement of grains, the newer high- 
speed conveying systems now handle chemicals, coal and 
a wide variety of industrial materials as well. Recent 
elaborations in these transport techniques have led to the 
development of dense-phase conveying (where, for ex- 
ample, 200 Ib. of finely powdered coal can be transported 
in the presence of only one Ib. of conveying air) and 
pneumatic stowing, used by modern mining technology to 
fill underground voids w ith crushed stones of relativ ely 
large size. 

Chemical engineers were probably among the first to 
realize that the gas-phase component in this binary system 
could be made to play a more active role than that of a 
mere conveying medium. They were instrumental, during 
the past twenty years, in the dev elopment of new pro- 
cesses in which the unique contacting properties of 
solids-gas flow could be made to effect either simple 
physical operations, such as evaporation and drying (as 
in flash drying) or chemical reactions, as illustrated by 
moving-bed or fluidized-bed applications too numerous to 
mention, 

It was inevitable that the spectacular achievements 
already obtained would cause a further intensification of 
research activities. As a result, two new contacting 
methods—the spouted-bed technique and the 
suspension technique—have recently been reported, and 


atomized 


new applications of solids-gas contacting are being de- 
veloped at an ever-increasing rate 

The advantages resulting from the use of finely- 
divided solid reactants or catalysts in chemical processing 
are many. Chief among them are the handling, 
the provision of large areas for and heat 
transfer, and the near-absence of internal diffusional 
resistance during reaction. The pulverization of the bulk 
raw material to yield particles small enough for efficient 


ease of 


surface mass 


contacting presented until recently a number of technical 
or economic problems, but remarkable improvements in 
comminuting machinery which have occurred during the 
past few years now make it possible to subdiv ide effec- 
tively and economically not only very hard minerals such 
as taconites and 
as wet bark, sawdust, peat and lignite. 


but even fibrous materials such 
Lhese results have 
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been achieved both by refinement of old designs (rod, ball 
or hammer mills) and by dev elopment of new principles 
of attrition (Aerofall mill, fluid erosion, explosive shat- 
tering ). 

As is commonly the case in many of the fast-develop- 
ing areas of chemical engineering, practical know-how 
and purely technical consi iderations have far outdistanced 
basic studies of the principles underlyi ing solids-g as con- 
tacting oper: ations. Even in the field of solids- -gas flow— 
which only involves considerations of momentum transfer 

the approach is still largely empirical and rests heavily 
on arbitrary friction factors, specific to each application, 
from which generalization 1s utterly impossible. It is 
perhaps surprising to report that a large number of work- 
ers, currently active in the field, are still following this 
line of attack. 

The discussion which follows was originally planned 
as a literature survey to record the results of experimental 
investigations of the conveying of solid particles by gas 
streams. Since in its broad aspects pneumatic conveying 
is a complex problem in fluid mechanics, the solution of 
which involves concepts of turbulence, flow fields, sto- 
chastic phenomena and boundary layer behavior, the 
search for experimental evidence had to be broadened 
beyond the usual limits of the chemical engineering 
literature, to include contributions from many diverse 
fields which include aerodynamics, fuel combustion, 
meteorology, acoustics, mining technology, industrial hy- 
giene and toxicology, geophysics and applied mathematics. 
That this search was highly productive is attested by the 
large number of pertinent references—some obtained with 
difficulty which have been consulted. The list is by no 
means complete. 

As the review progressed, it was noted that a good 
many of the experimental investigations reported had 
been aware of only a few of the many possible system 
variables, so that they dealt with a comparatively narrow 
range of conditions. As a result, most of the published 
data on solids-gas systems were restricted in scope and 
frequently conflicting. Also, the scarcity of fundamental 
information and the lack of uniformity of presentation 
made a comparative study of the experimental evidence 
extremely arduous. It soon became evident that a mere 
recital of a large number of apparently unrelated facts 
would simply compound the confusion. 

Although one of the purposes of the present work 
stil to record the published information as completely 
as possible, it was decided to place the main emphasis on 
critica] analysis and interpretation of the experimental 
evidence, using the microscopic in an attempt to unravel 
the various macroscopic phenomena involved. In view of 
the complexity of the inter-relationships existing between 
the numerous variables, this was no easy task. The sym- 
bolism and nomenclature in themselves presented quite a 
problem; an attempt was made at unification, which may 
not have been successful. 


Jo facilitate the presentation, the work has been 
divided into a number of self-contained articles, each with 
its OWN Nomenclature and bibliography, starting W ith the 
behavior of a single spherical particle in infinite, tur- 
bulence-free surroundings, and then going on to conside! 
all the various complicating factors, such as turbulent 
fields, effects of particle shape, roughness, acceleration and 
rotation, finally discussing the various aspects of multi 
particle systems. Thus the pertinency of some of the 
definitions and topics discussed at the beginning to the 
problems of solids-gas How may not become evident until 


much later. 
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In spite of its limitations, it is hoped that the present 
work will at least show the bolder outlines of this 
fascinating field and thus provide a better understanding 
of some of its fundamental aspects. 


FLUID PROPERTIES 


The two broad classifications of fluid flow—laminar 
and turbulent—have been extensively discussed in the 
literature and excellent texts are available giving full 
treatment to the more recent work in the field (* 4 5, 6, 7), 
A brief and at times oversimplified statement of some 
elementary concepts together with some less evident 
information which will be referred to in subsequent 
discussions will now be given. 


Laminar Flow 


In this flow regime (also referred to as viscous or 
streamline flow) the fluid moves in layers or lamina. It 
possesses a definite “viscosity” which is made manifest by 
the internal frictional resistance which it offers to phy sical 
deformation. Fluid movement in the vicinity of a bound- 
ary will cause the layers to slip over one another, and they 
will produce a shearing stress given by the conventional 
Newtonian relationship. The generalized movement of an 
element of fluid resulting from its reaction to the viscous 
and pressure forces which it experiences is concisely 
stated by the vector form of the Navier-Stokes equation 
for ce flow: 


Du/D6é@ = — (1/p) grad p + v \72u ay 


a Du/D6@ is the total derivative of the velocity vector, 
1, with respect to time, 6, and represents the total accele- 
ration of a fluid element due to changes of velocity in 

both space and time. F represents any body force (eg. 

gravity); this term may be neglected in the case of sub- 

merged flow with no free surface, and, in this case, the x 

component of Equation (1) may be written out in full 

as follows: 

ou bu ou du 1 6p 


ae Spee mp meee ape Tat 
60 6x by Oz p bx 


v 


Ou = 6’u_——s 2 
+ 7 + (2) 


6x? 62? 


From the law of conservation of mass, another im- 
portant equation—the continuity equation—can be derived: 


div(pu) = 6p; 60 (3) 


For an incompressible fluid, this expression may be 
written out in full: 


6u/dx + 6v/dy + dbw/éz = 0. (4) 


It should be noted that the assumption of incompres- 
sibility is valid for the gas phase in the immediate vicinity 
of a solid particle in most solids- gas systems. 


The Navier-Stokes equations are at the core of any 
theoretical description of a real fluid system, so that their 
specific solution would predict all of the desired flow 
properties for that system. The equations, however, 
usually contain more unknowns than the number of 
known pertinent equations for the system and their solu- 
tion presents almost insuperable mathematical complexity 
due to the non-linear terms. 


Simplifying assumptions and secondary behavioral laws 
have been successfully employed to obtain solutions for a 
limited number of cases. Advances in the mathematical 
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sciences coupled with the evolution of physical laws 
resulting from experiment may be expected to increase 
the number of specific fluid systems for which a solution 
is available. 

it is seldom justifiable in chemical engineering prob- 
lems to assume that the fluid is inviscid, that is to say, 
that it possesses zero viscosity. In certain cases, however, 
where the movement of the main mass of fluid, away 
from solid boundaries, is of particular interest, such an 
assumption greatly simplifies the equations of motion. If 
the viscosity terms are neglected in (1) and (2), the 
resulting expressions are called the Euler equations of 
motion. A detailed discussion of the solution of these 
equations—the subject matter of classical hydrodynamics— 
does not properly belong in this report. Suffice it to say 
that to integrate these equations, it is convenient to intro- 
duce a new function ¢, which is called the velocity 
potential, the derivation of which in any direction gives 
the velocity in that direction, or: 


ataname dale (5) 


=u; — =V; 


6x by 


If the flow is steady (8u/8@ = 0), the equations of 
motion can be solved directly along a streamiline, i.e., the 
line drawn tangent to the velocity vector at every point 
in the flow field. 


If the flow is steady and also irrotational, that is: 


dv du dus bw dw bv 
ab cee We tee Ree wey. ee oe ee 
bx dy 6z 6x dy & 


integration of the Euler equations for such a flow field 
(w hich is often called potential flow) yields the familiar 
Bernouilli equation (if gravity is the only external force 
acting on the fluid): 


(1/2)(V2/g.) + gZ/g. + p/p = constant ... (7) 


The two-dimensional analysis of many problems in- 
volving inviscid fluids may be represented by a network 
of orthogonal lines consisting of streamlines and of equi- 
potential lines. The latter are curves of equal velocity 
potential 4, while the streamlines are curves of equal 
stream function y. The latter is mathematically defined, 
for two dimensional incompressible flow, by: 


6y/dy = uand 6y~/dx = Vv ..(8) 


It should be noted that although the velocity potential 
concept is valid only for potential (or irrotational) flow, 
that of the stream function is equally applicable to viscous 
systems and is frequently used to map out complex 
laminar flow systems, such as the viscous flow around an 
object submerged in a moving fluid. 

If the flow is no longer irrotational, the three equa- 
tions (6) will not be equal to zero, but will represent 
twice the value of the angular fluid velocity, or the com- 
ponents of vorticity. Thus: 


bv/dx — bu/dy 2w, = vorticity in a plane parallel to the z axis 


(9) 
and the terms on the left of the three equations (6) are 
the vorticity components. The concept of vorticity is 
very important in boundary theory or in the 
ap proach to turbulent motion. Mathematically, vorticity 
is a vector oer obtained by taking the curl of the 
flow velocity, o1 


layer 


vorticity = curlu = \f Xu (10) 
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Physically it is a measure of the rotating motion associated 
with an infinitesimally small volume sample taken from 
the fluid. The locus of the centres of the circulation are 
often referred to as the vortex filaments or threads. 


Transition to Turbulent Motion 


As the velocity of flow of a laminar network is 
increased, the stability of the laminar layers diminishes 
and at a sufficiently high velocity a confused eddying, 
indicating turbulent motion, begins to be observed. 


There is as yet no clear understanding of the mechan- 
ism governing the transition from instability to full 
turbulence, although energetic inquiries towards this end 
have recently been reported (8% 19,11), Hsu (2) has re- 
cently proposed a plausible conjectured model of transi- 
tion to turbulence which is useful at this point in the 
discussion. The instability of a laminar layer is seen to 
produce a vortex filament, which moves relative to the 
ambient fluid. Since there is a drag force on the filament 
due to this relative motion, it is caused to accelerate. The 
acceleration causes increased vortex instability so that it 
breaks down resulting i in the generation of smaller vortex 
filaments which are in turn¢accelerated and distorted by 
ambient vortex filaments or by a shear flow which may 
be present, and this cascade process will be continued so 
that after several cycles there is a random orientation 
and distribution and size spectrum of the vortex filaments 
resulting in a field of random vorticity referred to as 
turbulence. The flow may return to the laminar state if 
the necessary instability conditions or other factors spe- 
cifically required at each stage of the cascade are not 
present to propagate the chain reaction. 

In his classical treatise, Reynolds “'3) found that the 
transition from laminar to turbulent motion of the flow 
in a pipe duct could be conveniently identified using the 
Up 


grouping of variables as a criterion. This dimension- 


less group is the ratio of a typical inertia reaction per unit 


U2 U 
volume a to a typical viscous force 7 It is a similarity 
parameter, and the greater its value, the smaller is the 


role play ed by viscous forces in determining the pattern 
of motion. 


The Reynolds grouping designated by the symbol Re, 
has been employ ed to predict whether fluid motion in a 
pipe is laminar or turbulent. The lower critical value of 
Re for duct flow is defined as the one below which 
turbulence cannot exist, ie. the viscous effects are suffi- 
ciently great to damp out and destroy fluctuations caused 
by momentary disturbances (!4), 

Typical of recent re-evaluations of Reynolds’ data is 

1 paper by Lindgren (13) who shows evidence that the 
claeaienl experiments used to determine the lower critical 
value of Re in a conduit—usually given as 2100—have been 
misinterpreted and that there does not exist a critical 
value in this region. His experiments further indicate that 
(1) Re alone does not determine the state of motion of 
fluids in pipe; (2) the transition from laminar to turbulent 
motion continues at least to Re 3300; (3) a sudden 
increase in resistance to flow at certain values of Re does 
not signify an instantaneous onset of turbulence; (+) the 
severity of the turbulence with high entrance disturbances 
decreases with the distance trav elled; this has been ob- 
served with flows up to Re 1,900,000; (5) in defining 
the state of motion in smooth cylindrical tubes, Re should 
be replaced by an expression which is a function of 
(DUp) (D) 

x  § 
entrance. 


, Where L is the downstream distance from the 


A comprehensive discussion of this subject is 
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available in the literature review presented by M. H. 
Smith (16), 

The transition region at which the flow in a pipe 
becomes turbulent may be considerably altered by the 
presence of a solid phase, and this will be discussed more 
fully in a later section. 


Turbulent Flow 

It is convenient to regard the turbulent motion of a 
fluid as made up of two separate components: (1) a 
mean forward motion of the fluid which is steady with 
time; (2) a superimposed secondary randomly fluctuating 
motion of the fluid whose time-averaged value is every- 
where null. Thus, if u, v, w and p are the instantaneous 
velocity components and pressure, the following typical 
relationship is obtained: 


EE hi ncc vith 
where u’ is the amount of fluctuation in the component 
velocity and 7% is its time average value, defined by: 

6 


er 


The time interval 0 — 6 is made sufficiently large to 
obtain a statistical sample so that a second sampling over 
this interval will not alter the average value u. 
It follows that: 
6 
8 


’ 


u = 


ON ME oo cece awn kee (13) 


Qi— 


~ 


to) 


Ie 


while the root mean square value of the fluctuation is 
given by: 
a 1 ; 
ly? 
VO Oly 


Similar expressions can be written for the other velocity 
components and for the pressure. 


° 


Turbulent flow can be crudely visualized as being 
made up of eddies which vary continuously in size and 
translational speed. The nature of turbulent movement 
is very complex and statistical methods are therefore 
necessary for its analysis. 

The root mean square values of the velocity fluc- 
tuations are referred to as the intensities of the turbulence 
components. They characterize the velocity within the 
eddies—packets of fluid which are constantly being formed 
and destroyed, their size varying from the order of the 
dimension of the pipe to v anishingly small “micro eddies” 

Whereas in viscous flow the viscosity depends on the 
interaction of fluid molecules, fluid stresses being trans- 
mitted by inter-mixing and collision of neighbouring 
molecules, in turbulent motion there is still interchange 
of matter but now on a molar scale, and there is a diffu- 
sion in which fluid packets from low velocity regions are 
transferred to high velocity regions with an accompanying 
retardation to the Jatter. 

The shearing stresses are therefore made up of two 
components: (1) the molecular action in the mean motion, 
and (2) the turbulent action in the secondary flow. The 
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latter are made evident by the so-called Reynolds shear 


stresses, —pu'v’, —pv'w’ and —pu'w’—three of six terms 
which occur in the turbulent flow form of the Navier- 
Stokes equations when the velocity is expressed as the 
sum of a mean and of a fluctuating component. The 
— three terms are the Reynolds normal stresses, —pu’, 
—pv"™, and —pw” 

There was an  eclgiind tendency to draw an erroneous 
analogy between eddy motion and molecular motion 
which involves the movement of rigid particles, resulting 
in a number of crude concepts which did, nevertheless, 
prove useful in providing a framework for further theo- 
retical and experimental work. Some of these expressions, 
which will arise in discussions in solids-gas flow, are 
listed below. 


An eddy viscosity y is defined by the expression: 
Ee NI EY inion eae io 0 ccmictanee (16) 


In this expression 7 is the total shear stress and is the 
sum of the laminar and of the turbulent shear stresses. 
The kinematic eddy viscosity is given by 1/p. 

Prandtl") defined a “mischungsweg” or mixing 
length \' the eddy counterpart of the molecular free path, 
at the end of which the volume element would transfer 
all its momentum to the adjacent fluid, so that, instan- 
taneously: 


F i INO ere ey 1616 Giger Subse Mace AUR (17) 


This implies that the eddy remains intact as a discrete 
entity over the length y. “Tay lor 8) in an early work 
presented a vorticity transfer theory, which sought to 
remove some objections to the Prandtl concept, by con- 
sidering a length |” over which the volume element trans- 
ferred all its excess vorticity rather than momentum so 
that, instantaneously: 


ae PGA) on xh ine es KER neal (18) 


The physical models upon which the above theories are 
based have since been shown to be entirely unrealistic. 
The eddies in a turbulent field are not discrete, and it 
is now considered naive to associate any of the length 
scales employed with a true physical picture of the 
turbulence. 

Taylor (1% 2021) recognized that eddy diffusion in- 
volved a continuous change of the eddy mass and size, 
there being a constant decay and loss of eddies through 
viscous shear, in equilibrium with a constant generation 
by the instability of the flow. This seemingly chaotic 
system was therefore treated by statistical methods, and 
was characterized by the correlation of the components 
of the velocity fluctuations at different space and time 
intervals. The correlation between two infinitesimally 
displaced points is at a maximum (defined as unity ) since 
the volume elements that would engulf one would be 
sure to include the other. The correlation would usually 
be expected to decrease as the intervals increased. Taylor 
defined a correlating coefficient Ry by: 


u’; u’s 


fe = safes » 
where u’ 


, and u’, are the velocity fluctuations at the 
points 1 and 2 separated in a . eee transverse to the 
fluctuations by the distance y. Similar correlation coeffi- 
cients can be ‘defined for the dictions x and z, and also 
for v’ and w’. 
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A characteristic length scale Ly is associated with Ry 
through the relationship: 


Le = | Ray coeue Othe 66 ee boise s (20) 


This term is referred to as a Eulerian scale. 

Obviously there are many Eulerian correlations pos- 
sible corresponding to the different components of the 
velocity fluctuations and to the three directions, with no 
simple relationship existing between them, so that a given 
flow system is only partially characterized by a given 
Eulerian scale. Triple and quadruple correlations are also 
employed. 

Consideration of the correlation between the velocity 
components of a fluid element separated by two intervals 
of time leads to the difficultly obtained Lagrangian scale 
L,, defined by: 


where 
Uglg + Ae 
‘2 2 
Vu @ Vu 6+ Aé 
There is no general relationship between the Eulerian 


and the Lagrangian scale (22) but limited proportionalities 
exist in certain cases. The mean square eddy diffusion 


Ry = 





distance X® is given by Taylor’s relationship: 
a 
1 dX? : 
Wt ROMP: cv ccateslere as (22) 
2 db | 
° 
1 dX? ; a 
so that > ig can be regarded as the eddy diffusion co- 
- £ , 


efficient. Its measure is often employed to study the 
effect of entrained solids on the turbulence parameters of 
the gas phase of a solids-gas system (24 24 25,26), Other 
statistical precepts introduced by Taylor include that of 
the turbulence spectrum, by means of which is indicated 
the relative turbulent energies associated with the various 
frequencies of velocity fluctuations which go to make up 
the total field of random vorticity. 

An _ idealisation removing ‘almost insurmountable 
mathematical difficulties is the concept of homogeneous 
oe turbulence, the study of which was begun by 

Taylor in 1935. Isotropic turbulence infers that the 
mean values of functions of the flow variables are not 
affected by the translation, rotation and reflection of the 
axes of reference. As a consequence, the intensity com- 
ponents in such flow are equal, or 


12 


YF YF We oped wens le ness ete) 


The concept of homogeneous turbulence may be 
loosely defined by the condition that the intensity com- 
ponents should not be a function of position in space. 
In other words, although the three intensity components 
are not necessarily equi al to each other, they remain un- 
changed at any point in the turbulent field. 

It should be noted that, although homogeneous, the 
turbulence may be either isotropic or anisotropic. Iso- 
tropic turbulence appears to be statistically the more 
probable state, and it was thought at one time that any 
anisotropic turbulent motion rapidly approached the 
isotropic state, but recent experimental work has shown 
that the approach to isotropy is very slow “), 

Batchelor 4) points out that although the kinds of 
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turbulent motion encountered in aeronautics, hydraulics 
and chemical engineering are usually more complicated 
than homogeneous turbulence, the study of the latter has 
practical utility in that an understanding of this simpler 
case will lead to comprehension of at least some of the 
aspects of non-homogeneous turbulence. On the other 
hand, there have been instances where investigators of 
solid-fluid flow have erroneously assumed isotropic tur- 
bulence to exist in a region of their experimental systems, 
and have then applied relationships derived from this 
assumption to obtain quantitative data. 

Although the literature abounds with experimental 
probings of the intercorrelation of the turbulence pro- 
perties, there has been little progress made towards the 
presentation of a universal theory of turbulence necessary 
if a solution to the Navier-Stokes equations for turbulent 
systems are ever to be solved. A typical illustration of 
the immaturity of our present state of knowledge is 
found in a recent paper by Betchov ‘?7) who points out 
that attempts to enunciate the mysterious processes w hich 
determine and lead to the production of vorticity lead to 
systems of equations which are always less than the 
number of unknowns. 


Fully Developed Turbulent Pipe Flow 


The turbulence structure of fully-developed particle- 
free pipe flow has been carefully ‘studied by hot-wire 
methods (2: 29 30,31), the hot wire anemometer being a 
direction-sensitive instrument whose electrically-heated 
fine wire reacts to the instantaneous convective cooling 

caused by the fluctuating flow (32, 33, 28), 

Laufer showed the pipe flow field to be divided into 
three concentric regions: (1) the wall vicinity in which 
the production and diffusion of turbulence, and viscous 
action play equally important roles; (2) the central 
portion of the pipe cross-section in which energy diffusion 
was found to be the important phenomenon, (3) the inter- 
mediate region in which the energy obtained by diffusion 
was dominated by the local rate of change of turbulent 
energy production. Typical turbulence profiles of the 
axial intensity as a function of ratio of the distance from 
7 pipe center-line r to the pipe radius r, are given by 

Laufer in Figures 1 and 2, the latter r amplifyi ing the inter- 
esting intensity peaks occurring near the pipe wall. Figure 
3 shows the somewhat similar behavior of the tangential 
and radial turbulence components but they clearly indi- 

cate that the flow is homogeneous and isotropic only over 

the centre fifth of the duct. The intensities in these plots 
were taken as the ratio of the occurring fluctuation to 
the friction velocity Ur defined by: 


: (%) 
Ue Rie bk a kc aeee ceees (24) 
a) an 


Turbulence intensities for pipe flow are usually quoted 


Vu’? 
U 


the axis; stated in this manner, the intensities at the axis 
would be approximately 2.7%. If the local results are 
given as the Jocal intensities, ie. using the local mean 
velocities as the denominator, this w ould indicate local 
turbulence intensity values which are as high as 60% for 
the pipe wall region. 

Similar, though less detailed results were obtained by 
Hall (34), Newman (35) and Sandborn (*)), the latter show- 
ing that the longitudinal turbulence intensity at the pipe 
centre was given by the general expression: 


Vu’? 


as the ratio where U, is the mean fluid velocity at 


= 0.144 Re-®.'. .. vc (25) 





— 
4 

} 
+ 





o) 
bo} 
bs 
+ 


“te 
Figure 1—Intensity of turbulence y,/2,[), versus distance 
from pipe wall (Laufer). 


Mickelsen ‘*) found that the Eulerian and Lagrangian 
scales at the pipe centre in fully -developed pipe flow 
were roughly related by the expression: 


Le = 0.6L; ; ; : (26) 


for turbulence intensities (reported as y ,’2) of 1.8 to 14 
feet per second. 


. FLUID RESISTANCE TO UNIFORM 
MOTION OF SINGLE SPHERICAL PARTICLES 


Spherical Particle Behavior in a Laminar Fluid— 
The Drag Coefficient 

The laws of motion of a single spherical particle in 
an undisturbed fluid stream (i.e. one in which there is no 
secondary motion) have been the subject of a great many 
investigations performed during the past two centuries. 

Newton ‘#7) in 1719 presented experimental data for 
the retardatory air force on hollow spheres which he 
dropped from the dome of St. Paul’s ¢ -athedral, and he 
derived the following equation which seemed to predict 
his results theoretically: 


R = (7/2) pr,?V? (27) 


Equation (27) also roughly approximates some of the 
results of more recent experience, although the premises 
on which the derivation is based appear to have little 
or no validity. This must be stressed, since it is often 


| ‘fy | _ 
o o + x Oo 6 | | 
- oe } } 
ery = 
} y | 
> | | 
7 , 
vi 4) t 
/U, i | | | | 
Ke | | 
| +0 500,000 +] 


° 50,000 





1.6F-O ; } i { ' ~| 
“a | 
yor ot | | 
etd oS ag - | | | 
i, | — 
1 | 
4 Lacan” 
| 
‘cae. Pk 
Q \ 2 3 G 5 6 7 
wiiedlines 
i 


Figure 3—lIntensity of turbulence y¥\2 (| and y\/2/| 


7 
versus distance from pipe wall (Laufer). 
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Figure 2—Intensity of turbulence ¥ |)/2/[), versus distance 


from pipe wall in wall vicinity (Laufer). 


assumed to be a successful solution to the problem. It 
will be subsequently shown that it is no paradox that 
subsequent writers have been unable to propose a com- 
plete theoretical solution for the flow regime for which 
Newton’s law seems to be valid. Zahm (#8) gives the fol- 
lowing summary of Newton’s assumptions: (1) the fluid 
is a discontinuous rare medium consisting of infinitely 
small equal particles, freely disposed at equal and great 
distances from each other; (2) these fluid particles, in 
uniform rectilinear motion, strike only the front part of 
the solid body immersed in it, leavi ing a vacuum behind 
extending to infinity and bounded by parallel streamlines 
grazing the body; (3) the particles striking the hard 
smooth front retain their tangential velocity component 
and lose their normal component if inelastic, or have it 
reversed if elastic. In all cases before and after rebound, 
they continue in straight lines without mutual interfer- 
ence, ie. they have substantially infinite free paths. 

The laws of resistance derived from these premises are 
absolutely valid for such a system, but this system differs 
vastly from a natural fluid. The Newtonian conception 
of fluid resistance has been replaced by the hydrodynam- 
ical conception, according to which ‘the resistance con- 
sists of pressure differences and frictional stresses arising 
from the fluid around the body (39), 

The phenomenon of fluid resistance is extremely 
complex and no simple formula can express the relation- 
ship between the fluid flow and the force on the spherical 
particle. Limited success has been attained for restricted 
conditions; however, results covering a wide range of 
conditions can only be given graphically. 

A most important simplification was discovered by 
Lord Rayleigh ‘#) when he found that a single graph 
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Figure 4—Standard drag coefficient curve. 
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could be used to present data for the motion in a laminar 
fluid of spheres of all sizes, and fluids of all densities and 
viscosities. One of the coordinates uses a ratio of inertia 
reaction to viscous reaction on the sphere, and it is 
expressed as a modified Reynolds Number (d,Vp)/u-— 
where d, is now the diameter of the sphere. The other 
coordinate is a force coefficient R/(1/2)pV?A, in which 
R is the resistance of the fluid to the motion of the sphere 
and A, is the cross-sectional area of the particle. This 
dimensionless group is generally referred to as the co- 
efficient of drag and is designated by the symbol “C,”. 

The results of the many investigations were collected 
and graphically correlated using the drag coefficient and 
the particle Reynolds Numbers as coordinates. This has 
been presented by Castleman (41), Eisner (42), Muttray (4), 
Rouse (#4), and most recently in a plot by Lapple and 
Shepherd (#5). 

Plotted on logarithmic coordinates, as shown in Figure 
4, the drag coefficient falls linearly from a value of 
100,000 at a Reynolds Number of 0.00025 to a value of 
240 at Re, = 0.1. This is the Stokes Law region, the 
significance of which will be discussed further on. It then 
decreases more moderately to a minimum of 0.38 at Re, 
= 5000 and rises gently to a plateau of 0.47 at Re, = 
100,000. It falls off sharply to a value of approximately 
0.1 at Re, = 300,000 and again rises slowly with the 
Reynolds Number past that point. 

At first glance it would appear that the determination 
of the drag coefficient would require only a minimal 
program of experimentation, involving the efforts of only 
a handful of workers. However, the wide discrepancy in 
the original results emphasizes the difficulties in Ubtaining 
the required idealized environmental conditions and in 
making accurate velocity measurements. The authors have 
therefore deemed it pertinent and informative to list, in 
a chronological sequence, the more significant contribu- 
tions made towards the evolution of the standard resis- 
tance curve for spheres. This curve is meant to apply 
only to smooth non-rotating spheres, moving at a constant 
relative velocity in an infinite fluid free from any dis- 
turbances. Conditions of observation not complying with 
these restrictions usually result in values which deviate 
from the curve in some instances, as will be seen, by as 
much as several thousand percent. 


In 1867, Rittinger r (46) sought to verify Newton’s “Law’ 
by inv estigating particles whose Reynolds Numbers were 
in the vicinity of 10,000. He concluded that Newton’s 
formula was not valid since constant k, in the relationship: 


R = ko(w/2)r3,V?........ oa (28) 


was not unity, but varied with the particle Reynolds 
Number. 

A fuller Reynolds Number spectrum of from 0.009 to 
8200 was investigated i in 1900 by Allen (*), who developed 
a stereoscopic photography technique to measure the 
velocities of essentially spherical air bubbles and spheres 
of amber, paraffin and steel moving in water and aniline. 
His values of C,, are at times low, possibly because the 
terminal velocities may not have been attained. They are 
also complicated by boundary or wall effects. : 

Arnold in 1911 (4) reported the fall velocities of rose 
metal spheres of diameter 0.013 to 0.14 cm. in colza oil, 
corresponding to a particle Reynolds Number of 0.002 to 
2.4; however, his results are complicated by a wall effect. 


Since it is easier to measure the velocity of a fluid than 
that of the particle, a solution to the drag coefticient 


problem was sought by somehow measuring the force on 


spheres suspended a laminar fluid stream. The errors 
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due to unsteady relative motion which plague the fall 
velocity studies are eliminated if the fluid velocity is held 
constant. Other difficulties develop, however, due to the 
inherent instability of a laminar fluid moving in a large 
size duct, since its Reynolds Number would be in excess 
of 2000. 

Becker in 1907 (4°) measured the resistance of spheres 
suspended in an air stream, but his results are not valid 
because they are overshadowed by the effect of the 
supporting wires. There are two principal sources of 
error here: the evaluation of the fluid force on the sup- 
port mechanism is difficult—this is no simple matter, since 
the presence of the sphere may alter this force; secondly, 
the perforation by the supports of the complex flow field 
in the immediate vicinity of the sphere may in some cases 
engender unexpectedly large changes in the drag. This 
effect will be discussed more fully in a followi ing section. 
Similar difficulties were encountered by Constanzi who, 
in 1912 5°), reported data obtained for spheres suspended 
in a moving column of water. 

Meanwhile work with spheres moving in stagnant 
fluids was continuing with the help of improved velocity 
measuring techniques. Eiffel in 1912 5" obtained resis- 
tance coefficients for spheres dropped from the Eiffel 
Tower, over a path long enough to allow the attainment 
of a constant velocity. These values were later confirmed 
by Prandtl in an artificial current of air at the Goettingen 
wind tunnel (52), - 

Shakespear in 1914 (53) obtained very consistent results 
when measuring the fall velocity of weighted celluloid 
spheres in an air tower, by making use of an electrical 
timing system. His system was, however, incapable of 
investigating small Reynolds Numbers since it would 
have required extremely small and at the same time 
accurate spheres of low density. This region was there- 
fore almost exclusively studied with a liquid as the fluid. 


Schmidt in 1920 (54) recorded the velocities of wax 
spheres i in water, as well as the rising velocity of balloons 
in air. His results are complicated by acceleration effects 
which will be discussed more fully in another section. 
In 1921 Krey (55) conducted an extensive series of tests 
with lead, glass, brass, steel, clay, hard rubber and wooden 
spheres in water. These covered the Reynolds Number 
range from 100 to 30,000, which he extended to a value 
of 1 x 10° by measuring the force on a wooden sphere 
suspended in an air tunnel. 

The great surge of activity in the field of aeronautics 
which developed at the close of the first World War 
produced a need for air tunnels in which test models of 
aircraft would behave as they would in the effectually 
laminar conditions of the free atmosphere. The obstacles 
to be overcome were indeed the very same as those which 
caused error in the suspended-sphere experiments. It 
became evident that the turbulence present in the flow of 
the early wind tunnels caused the behavior of the model 
aircraft to be radically different from that of the proto- 
type in the atmosphere, so that the success of a modern 
wind tunnel design is in part measured by its ability to 
maintain streamline flow at high tunnel Reynolds Num- 
bers. The fluid movement is so adjusted that the tur- 
bulence which originates at the walls of the duct is not 
allowed to permeate the flow in the core (which contains 
the test section)—as would be the case in fully developed 
pipe flow. A wind tunnel which fulfilled these require- 
ments was constructed at Goettingen, and in 1922 
Wieselsberger % began a series of dr ag coefficient trials 
in which he suspended (with a minimum of error) alum- 
inium spheres of 1 to 28 cm. in diameter to get particle 
Reynolds Numbers of from 700 to 770,000. ‘His results 
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generally coincide with the modern standard curve, and 
this is a tribute to the uniformity and lack of turbulence 
in the Jaminar air stream produced at the Goettingen 
tunnel. 


Some rather novel methods were employed to extend 
the drag coefficient range. Richardson in 1923 (°*) measured 
the resistance of air to the movement of spheres by firing 
them upwards from the unrifled barrel of a cannon. His 
results are higher than the now accepted values because 
of acceleration effects. Bacon and Reid in 1924 (58) ob- 
tained results for very high particle Reynolds Numbers 
(2.5 x 105 to 1.5 x 10°) by observi ing the motion of wood 
and brass spheres dropped from a plane. 


Very accurate results for low Reynolds Numbers were 
obtained by Liebster and Schiller (1924) using photo- 
graphic methods incorporating a timing disc (59), These 
experiments were performed with lead and steel spheres 
in liquids of various viscosities, and they covered the 
Reynolds Number range of 0.12 to 2000. 


Knodel in 1926 ‘°) used a pendulum method to meas- 
ure the resistance of aluminum spheres of radius 0.146 to 
0.401 cm. suspended in a laminar stream in a relatively 
narrow pipe. He developed an ingenious procedure to 
compensate for the complications introduced by the 
presence of the suspending mechanism. His results’ agree 
with the resistance values given by Liebster and Schiller. 

One of the most important contributions towards the 
evolution of the standard curve was submitted by Lunnon 
in 1926 who performed a lengthy series of trials with 
spheres of various densities and diameters dropped in air 
(61) and in water ‘®), and he was able to reproduce the 
results of Wieselsberger’s air tunnel tests. This proved 
effectively that the fluid resistance encountered by a 
spherical particle depended only on the relative move- 
ment of the solid with respect to the fluid, and not on 
the absolute values. Lunnon also investigated the effect of 
the acceleration of the particles on the drag coefficient, 
and this will be discussed at length in another section. 

Air tunnel drag tests done in the Langley Memorial 
Aeronautical Laboratory in 1929 by Jacobs ‘®%) confirmed 
that specific drag coefficient value is obtained at a given 
particle Reynolds Number, irrespective of what combina- 
tion of variables, ice. kinematic viscosity, velocity and 
sphere size, is used to obtain that value of the Rey nolds 
Number. 

There are more recent works such as that by Moller 
4) in which measurements generally confirming the drag 
coefficient curve have been made as a secondary issue. 
Minor corrections may result from the work of Wiselius 
65) who in 1947 found minimum drag values of 0.05 for 
the Re, 3 x 106 region. 

The validity of the standard curve for the sensitive 
region of 10 < Re, < 100 has been questioned recently by 
MicKnowa ¢ 6). who cites experimental evidence which 
indicates that the results in this region are still suffering 
from wall effects. 

The studies just discussed were primarily directed 
towards the measurement of drag coefficients for 
idealized conditions. The large amount of experimentation 
studying the—at times—over-riding effects of fluid tur- 
bulence and boundary conditions, as well as particle 
acceleration, rotation, surface roughness and shape, will 
be presented after a more detailed discussion of the aero- 
dynamics of the idealized case. 


Approximate Equations for the Drag Coefficient 
I-mpirical expressions of limited accur acy giving either 
the drag coefficient or the dr: ag force on a sphere for 
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restricted Reynolds Number regions are available, the 
following of which are representative. 

Schiller and Nauman ‘*?) found the following relation- 
ship to hold for Re, < 700: 


2 


< 


1 
Cp = 2 CL = OS ‘Rey 0)... cee ce ccG9) 


Re, 


Kliachko (65) used the approximate relationship: 


24 ; 
Cp = dG ik cvs ecanecas (30) 
Re, 
Langmuir and Blodgett ‘*®) employed: 
Re, ms : : 
Cp. 24 = 1 + 0.197 Re,°-® + 0.0026 Re,'-38...... (31) 


with good accuracy for the range 1 < Re, < 100. 

Budryk (70) felt that in the ‘ ‘intermediate range” of 
particle movement, (i.e., in the Reynolds Number range 
from 0.5 to 2,000) the drag force should be made up of 
the sum of the resistances given by Stokes Law and 
Newton’s “Law”, giving the relationship as: 


R = 6wur,V + ks (=~ ) r,%pV?..........06- (32) 


where k, is a constant associated with Newton’s Law. 
This pseudo- -theoretical equation deviates considerably 
from observed values. On the other hand, Froessling ‘ 71) 
and later Gumz ‘7: 73) successfully applied the following 
empirical relationship for this region: 

sedans Shee (33) 

RR TRC e 3: 
where Rg, is the resistance given by Stokes Law and Ry 
is a residual drag force found by using the following 
pseudo-drag values: 


C’p = 0.08741 for 0 < Re, < 8 
C'p = 0.3 for 8 < Re, < 300 
C’p = 0.28 for 300 < Re, < 2000 


Flow Regime Variation with Reynolds Number 

A brief description will now be given of the alterations 
which occur in the flow field about an immersed sphere as 
the particle Reynolds Number increases, in order to define 
the various flow regimes which will be subsequently dis- 
cussed in detail. 

The velocity of the fluid moving past the sphere “= 
be zero along ‘the surface of contact with the solid. 
very low Reynolds Numbers, the velocity will increase 
gr radually as the distance from the surface increases, to 
attain eventually a value which is equal to the stream 
velocity. This flow has up- and down-stream symmetry. 
The elements of the fluid meeting the front of the sphere 
are slowly accelerated sideways, and the inertia effects 
are too small to cause a time lag in the closing up of the 
flow behind the body. This visualization corresponds to 
the Stokes Law region which ends at a Reynolds Number 
of approximately 0.1, ushering in a “transition” region 
which has customarily been assumed to extend up to a 
Reynolds Number of 20,000. 

A separation process in the boundray layer is first 
hinted at a Reynolds Number of approximately one, at 
which the streamline patterns associated with Stokes’ 
regime become deformed and curl up to form a stationary 
vortex ring imbedded in the boundary layer at the rear 
of the sphere (4), ‘This ring resembles a smoke ring and 
it has strict axial symmetry. The velocity of circulation 
of the vortex fluid in contact with the rear portion of 
the sphere appears to be much slower than that of the 
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fluid in contact with its forward portion (75), The ring 
grows in size as the Reynolds Number is increased, the 
vorticity generated in the ring all the while being balanced 
by that which diffuses downstream from it. 

' As the Reynolds Number increases in the transition 
range, the velocity gradient in the boundary layer be- 
comes more severe. Since the fluid elements meeting the 
front of the sphere are accelerated sideways, this increase 
in the fluid velocity is accompanied by a decrease in 
static pressure, i.e. a negative pressure gradient, around 
the surface of the sphere towards the equator and in the 
direction of flow. The velocity of the fluid just outside 
the boundary layer begins to decrease beyond the equator, 
and the decrease in velocity now produces a positive 
pressure gradient. In that region, the motion of a thin 
stratum of fluid adjacent to the wall of the sphere and 
wholly inside the boundary layer will be influenced by 
three forces: it will be pulled downstream by the viscous 
action exerted by the stream above it; it will be retarded 
by friction at the sphere surface; and it is further retarded 
by the positive, or “adverse” pressure gradient. This pres- 
sure gradient increases with an increase in the sphere 
Reynolds Number and the energy and momentum of the 
fluid in the boundray layer become insufficient for it to 
force its way for very long against it. The fluid in the 
boundary lay er is therefore eventually brought to rest, 
with respect to the sphere surface, and farther on next 
to the wall. a slow back-flow in the direction of the 
pressure gradient may set in. The forward boundary 
layer stream will then leave the sphere surface, and the 
latitude at which this occurs is referred to as the point of 
separation. Figure 5 shows a tracer photograph of the 
flow in this region. 

At a Reynolds Number of about 500 (4), for a stream 
free from turbulence, a point of imbalance is reached 
and the vorticity diffusion rate downstream can no longer 
keep up with that generated in the vortex ring. At this 
point, referred to as the lower critical Reynolds Number, 
the vortex system breaks away from the body and feeds 
into a flow regime at the rear of the sphere which i 
bounded by a free shear layer, i.e. the separated beendury 
layer, and is referred to as the wake. The instability at 
the lower critical Reynolds Number does not correspond 
to that occurring at laminar-turbulent transition. As the 
Reynolds Number increases, there is a laminar-turbulent 


transition in the free shear layer which draws nearer to 
the separation point, the position depending on the free 
stream turbulence (7°), 

The transition region is traditionally assumed to ex- 
tend to approximately 20,000, after which the wake is 
supposed to be made up of completely random flow systems 
(77), Recent investigations point out, however, that there 
is always a periodicity associated with the turbulence in 
the wake, with the proportion of turbulent fluctuations 
gradually taking the upper hand as the Reynolds number 
increases (78), 

In the vicinity of Re, = 200,000 the forward part of 
the boundary lay er attached to the body becomes unstable 
and enters into a transitional state tending towards tur- 
bulence, producing an unsymmetrical turbulent wake in 
an unstable and shifting position. At the higher critical 
value of the Reynolds Number (Re, = 300,000), the 
separation circle moves suddenly to the rear, and there is 
a sharp drop of C, from 0.40 to 0.1. At still higher values 
of the Reynolds Number, the boundary layer flow be- 
comes turbulent and the wake, now relatively narrow, is 
also turbulent and statistically stable. 

The resistance of a sphere to motion in a fluid is 
spoken of as being the sum of a frictional effect, or skin 
drag, and a pressure gradient effect referred to as the 
form drag which is intimately associated with the mech- 
anism of flow in the wake region. The relative contribu- 
tion from each type of resistance varies with the flow 
regime, the former accounting for two-thirds of the total 
drag i in the Stokes Law region and its contribution steadily 
diminishing toa negligible \ ralue as the particle Rey nolds 
Number is increased beyond this range. 

The ability to describe experimentally and theoretic- 
ally the fluid ‘behavior about a spherical particle seems to 
decrease with increasing Reynolds Numbers. Successful 
solutions to the Navier-Stokes equations have been giv en 
only for very low values and the degree of approximation 
increases as the lower critical Reynolds number is ap- 
proached. The transition region has not been sufficiently 
probed by experiment, and has received only partial semi- 
theoretical treatment by boundary-layer theory methods. 
The developed wake regime has been even less fully 
described, and because of its complexity theoretical treat- 
ments deal only with the inter-relationships between 
arbitrary parameters. 





Figure 5—Flow around submerged body in the vicinity of the lower critical Reynolds number. 
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Our detailed discussion will follow this pattern, start- 
ing with the successful mathematical solutions available 
for the lower Reynolds Number region and ending with 
the most recently proposed concepts of wake flow 
behavior. 


ltl. THEORETICAL TREATMENT OF FLUID 
RESISTANCE TO SPHERE MOTION 


Stokes’ Solution 

A theoretical analysis of the problem of sphere motion 
should provide solutions for two major considerations: 
the first would permit calculation of the drag force for 
all Reynolds Number values. The second would provide 
a complete description of the flow field, so that the fine 
details of the changes occurring with increased velocities 

which can often not be discerned by experiment—w ould 
be revealed. In turn, the validity of a presented solution 
could be evaluated by a comparison with the available 
experimental evidence—the dr ag criterion being the less 
sensitive one, and in some instances being unreliable as 
was illustrated by the Newton’s “Law” case. Successful 
solutions would ultimately be extended to the less ideal- 
ized systems having boundary conditions imposed by 
shape, surface roughness, rotation, free stream turbulence, 
etc. As will be seen from what is to follow, incorporation 
of these additional complexities in a theoretical treatment 
will not be possible in the immediate future. 

In 1850 Stokes ‘7% succeeded in solving the equations 
for the resistance encountered by a sphere moving very 
slowly in a fluid, and this is the only case where a full 
theoretical solution has been obtained. The troublesome 
inertia terms of the Navier-Stokes equations may now 
be neglected so that it is reduced to the readily solvable 
form: 


grad p = v /*u (34) 
The net force on the sphere is given by 
Rs 3rdpyl (35) 


two-thirds of which are shown to be due to the shearing 
stress exerted by the fluid on the sphere surface, while 
the remaining third is due to pressure distribution. The 
resulting drag coefficient 24/Re corresponds exactly to 
experience up to a Reynolds Number of 0.1. 

This solution, considered a triumph of applied mathe- 
matics, is given in a simpler vector analysis form by 
Page and represents the only successful fully theo- 
retical derivation for the drag force on a particle im- 
bedded in a fluid to date. Stokes’ analysis provides a good 
aesc mpuon for the ambient flow exc cept in the region away 
from the sphere where the fluid movement approaches 
that of the free stream 


Oseen’s Approximation 


Neglect of the inertia terms becomes increasingly 
unportant at Reynolds Numbers above 0.1 and the Stokes 
¢ becomes correspondingly deficient. Whitehead 
attciipted unsuccessfully to calculate the flow pattern 
for hnite Keynolds alues by assulliing an expansion of 
the flow in terms of Re. Oseen pointed out that the 
result could not satisfy the boundary conditions, and he 
en proposed an approxunation to the Naviei Stokes 
ich onsidered the mil ing inertia terimis to 
) portant on} In that region of the field approa¢ hing 
iform flow. For this region the Navier-Stokes equations 
edaguced to the near form 
138 


1 ‘ 
U.gradu = —-p+p Vu........ ee 
p 


where U is the free stream velocity vector. The Oseen 
equation has often been given in a different form resulting 
from representation of the velocity distribution as the 
sum of a constant velocity U and a disturbance velocity 
u ‘), but Proudman and Pearson (87) point out that al- 
though it is formally the same as Equation (36) it is 
conceptually incorrect and has led to misleading and 
erroneous conclusions. They also feel that although 
Oseen’s method of solution had seemed to give drag 
results which correspond to experiment, it could not be 
considered superior to that of Stokes for the derivation 
of the drag. 

Goldstein 83) in 1929 succeeded in providing a com- 
plete solution for the Oseen approximation, and the 
following drag coefficient expression resulted: 


Cp = (24/Re,) (1 + 3Re,/16 — 19Re,2/1280 + 
71Re,3/20480-. ... .(37) 


This yields values which are within 5% of adi at 
Re, = 1 and 20% in excess at Re, = 10. Even more 
important are the flow patterns w hich are predicted by 
the solution, since the sphere now enters a Rey nolds 
Number region which sees the emergence of w ake phe- 
nomena together with the dev elopment of the complex 
wake structure. It was not until 1950 that the massive 
computations involved were performed by Tomotika and 
Aoi ‘), Their calculated pattern of streamlines clearly 
mdicaved the formation of a fixed vortex ring at the rear 
of the sphere at Reynolds Numbers as low as 0.1. In 1955, 
however, Pearcy and McHugh (*4) armed with a digital 
computer repeated these calculations, and were surprised 
to find that the solution did not show the faintest sugges- 
tion of the formation of the vortex ring even at Rey nolds 
Numbers as high as 10, where these rings had already been 
observed experimentally. Aside from pointing out serious 
errors in the work of Tomotika and Aoi, the calculations 
showed that the Oseen approximation indicates the forma- 
tion of a strongly marked wake at Reynolds Numbers as 
low as one, and a boundary layer type of flow appeared 
on the sphere surface in the vicinity of Re, = 10. Pearcy 
and McHugh’s results began to raise serious doubts as to 
the validity of the Oseen approximation for the flow about 
a — since it seemed that it was not c: apable of indi- 
cating a vortex no matter how high the Reynolds para- 
meter. 


Stewartson (55) felt that since the Oseen equation 
would have to be used as a guide in the eventual solution 
of the exact Navier-Stokes equations, he would investi- 
gate its validity by inspecting the flow which it would 
indicate at high Reynolds Numbers. The solution maps 
out flow lines indicating a boundary layer on the forward 
portion, and a wake extending to infinity at the rear in 
the form of a cylinder with its axis par: allel to the flow. 
The equations indicate that the wake is bounded by a 
shear layer, the flow inside and at infinity moving with 
the free stream velocity. At no time is a vortex indicated. 
This of course does not correspond to the actually ob- 
served system, The Oseen approximation of the Navier- 
Stokes equations is unacceptable as a description of the 
flow at high Reynolds Numbers since its solution indicates: 
(a) a bounds ary layer at the sphere front of order v, i.e. 
0 (v) whereas it is observed to be 0 (v%)—as indicated 
by the Navier-Stokes equation, (b) the Oseen-derived 
boundary layer does not break down when the skin 
friction vanishes on the body surface whereas in 1948 
Goldstein '* obtained a Navier-Stokes solution agreeing 
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with the observed fact that it does end and leave the 
surface of the body when this occurs; and (c) the solution 
for the pressure drag coefficient for Up/2u — © gave a 
value of 0.35. This, coupled with the Tomotika and Aoi 
solution for the Oseen approximation indicating that for 
all values of Up/2u the skin friction is § of the total drag, 
would indicate an obviously incorrect drag coefficient of 
1.05 as Up/2u (and therefore the Reynolds group) ap- 
proach infinity. 

In the light of the above, the Oseen approximation 
cannot, as it stands, be employed to further the theoretical 
advance towards a description of the flow for sphere 
Reynolds Numbers encompassing the formation and 
evolution of the wake regime. Proudman and Pearson (8?) 
have very recently presented an interesting solution for 
the stream functions for small particle Rey nolds Numbers 
by obtaining separate expansions of the stream functions 
for the region yen to and far away from the body 
immersed in the flow, referred to as the Stokes and Oseen 
expansions respectively. For the first expansion, the 
boundary condition dictates that there must be no slip 
at the surface, while for the second the stream must be 
uniform at the boundary. Although unique solutions could 
not be obtained, the authors succeeded in using a “match- 
ing” procedure which produced the alternate terms of 
each expansion. The solution was applied to the flow 
around spheres and cylinders, and for the former the 
following expression for the drag coefficient was obtained: 


Cp = (24/Re,) [1+ (3/16)Re, + (9/160)Re,? log Re, + 
O(Re,?/4)] . (38) 


The flow patterns about the sphere resulting from the 
solution have not as yet been calculated, so it is not known 
whether they would predict the vortex formation. 


Kawaguti’s Approach 

An — attempt to obtain an approximate solu- 
tion of the Navier-Stokes equations for Reynolds Num- 
bers in the vicinity of the lower critical value has recently 
been reported by Kawaguti ‘**: 5"), He employ ed a Galer- 
kin method which involved the assumption of various 
functional forms with unknown parameters for the stream 
function y, the unknowns being solved for by the bound- 
ary conditions and some simplification of the Navier 
Stokes equations. Many forms were attempted and the 
following form was successful for the range 10 < Re 
SO: 


v (1/2r? + Aj/r + Ao/r? + Az/r* + A,/r*) sin®?é 4 
(B,/r + Bo/r? + Bs/r? + By/r*) sin’@ cos 6 (39) 


where r and @ are the polar coordinates and 


120 72 
\, { \ i 
29 29 


95 17 
Ag 1 Aa is 
58 29 
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B, = — 9 Bi; 
5 

By =a 9 B; 


Values for A, and B, were obtained and the solution gave 
drag coefficients which are in excellent agreement with 
the experiment values in this range. Kawaguti reasoned 
that if Equation (39) is indeed correct, then the addition 
of a mathematically expressed small disturbance or per- 
turbation of the proper form to the equations would yield 
the value of the lower critical Reynolds Number. The 
disturbance would die out for any Reynolds Number 
smaller than this value and would increase for larger ones. 
In this manner a lower critical Reynolds Number of 51 
which Kawaguti compared with an experimental Reynolds 
Number of 100 was found. This comparison is approximate 
since the lower critical Reynolds Number occurs in the 
vicinity of 500 for a turbulence-free stream; all in all, 
however, this creditable approach shows promise. 

Kawaguti’s solution is not completely theoretical, since 
it relies on a comparison with experimental data for the 
selection of a proper mathematical form for the stream 
function. The work points out, however, that the exact 
solution of the Navier-Stokes equations together with 
the knowledge of the proper perturbation function could 
be expected to give the lower critical Reynolds Number 
as a function of the free stream conditions. 


Numerical Solutions 


Approximate numerical methods using finite ditference 
approximations have been employed in attem pts to obtain 
solutions at specific Reynolds Numbers ‘%: *); however, 
in order to obtain workable boundar vy conditions the 
sphere is usually assumed to be in a flow which is bounded 
by a wall w hich invariably exerts a marked stabilizing 
influence on the sphere velocity field. The most recent 
example of this approach is found in a paper by Jenson ‘#) 
in which the ambient fluid is contained in a evlinder 
having six times the diameter as the contained sphere, and 


the solutions which were obtained for Revnolds Numbers 


of 5, 10, 20, and 40, show separation angles which are 
consistently lower than the experimental results reported 
by Garner %), and wake dimensions which are smaller 
than those observed by Taneda ‘°°. 

In summary, although much effort has been expended 
in the study of the action on an immersed sphere by a 
turbulence-tree flow, a theoretical description of the 
important region above a Reynolds Number of unity, 
where the inertia terms of the Navier-Stokes equation 
has vet to be attained. The 
recent increased activity in the 


begin to exert an influence, 
field has been coincident 
with the evolution of tmproved mathematical methods 
and the application of computing devices to cope with the 
ponderous numerical calculations which are required to 
plot the tlow  tields 

\ subsequent paper will co uplete the discussion of the 
system of a smooth sphere moving in steady rectilinear 
motion through an essentially infinite and undisturbed 
thuid, by dealing at length with the many theoretical and 
experimental probings of the very unportant wake region 

Phe idealized boundary conditions thus far reviewed 
are seldom realized wo industrial solids-gas systems, and 


the at times over-riding ettects of non-steady motion, 


rotation, shape, roughness, wall proxunity and high rela- 


} 


tive intensities of turbulence will be discussed 
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Nomenclature 


Any consistent set of units may be employed. Those listed are 
merely illustrative 


Ap = characteristic particle cross-sectional area, ft.? 
Cp = particle drag coefficient, dimensionless. 
Cr = pseudo-drag coefficient used by Gumz, dimension- 
D J £ ) 
less Eqn. (33 
D = duct diameter, ft. 
d, = characteristic particle diameter, ft. 
F = body force, lb. force. 
i a a ts f oravity. ft./sec.2 
y = local acceieration O! gravity, It./sec, 
g. = conversion factor 32.17 (lb. mass) (ft.)/(Ib. force) 
(sec.?). 
ko, ks = Newton's ‘“‘Law”’ correction coefficients, dimension- 


less. 


downstream length, ft. 


L; = Eulerian length scale, ft. 

Ly = Lagrangian length scale, ft. 

I’ = Prandtl mixing length, ft. 

‘ig = Taylor mixing length, ft. 

R = particle resistance, lb. force. 

Ry = residual particle resistance given by Gumz, lb. force, 
Egan. (33). 

Rs: = particle resistance given by Stokes law, Ib. force, 
Egn. (35). 

Re = system Reynolds Number, dimensionless. 

Re, = particle Reynolds Number, dimensionless. 

Ry = Eulerian correlation coefficient, dimensionless. 

Re = Lagrangian correlation coefficient, dimensionless. 

i = distance from duct center, ft. 

r = pipe radius, ft. 

> = particle radius, ft 

U = mean stream velocity, ft./sec. 


uniform stream velocity vector. 


= velocity at duct axis, ft./sec 





U, = friction velocity, ft./sec., Eqn. (24) 
u = velocity vector, ft./sec. 
u, V, W = velocity components of U in Cartesian coordinates 
x, y and z, respectively, ft./sec. 
u’, v’, w’ = instantaneous components of velocity fluctuation in 
x. ¥, £, 3t./ sec, 
u = time average value of u, ft./sec, 
V = bulk fluid velocity, ft./sec. 
xX? = mean square eddy diffusion distance, ft.? 
= height, ft 
6 = partial differential operator 
7 = eddy ViISCOSILYy, Ib. mass/ft. sec 
6 time, sec. 
K viscosity, ib. Mass/It. sec 
v kinematic viscosity, ft.*/sec 
p fluid density, lb. mass/ft.* 
T shear stress, lb, force/t{t.? 
gd velocity potential, ft.2/sec., Eqn. (5) 
Y = stream function, {t.*/sec., Eqn. (8) 
Ww rotational velocity, radians/sec 
vector Operator del 
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Thermodynamic Properties of Some 
Organic Compounds 


I. PVT Relationships and Calculated Thermodynamic Properties for Normal Butanol* 


L. W. SHEMILT?, R. W. 


Pressure-volume-temperature measurements on 
n-butanol have been carried out using the Kay type 
of pressure-block apparatus and glass capillary cell. 
Determinations were made isothermally in the satu- 
rated vapor and liquid and superheated vapor regions 
to above the critical point. Compressibility factors 
have been calculated. Constants in the Beattie-Bridg- 
man equation of state have been determined from 
250° to 300°C. at pressures up to 75 atmospheres. 


The thermodynamic properties of normal butanol 
have been calculated over a temperature range of 
0°C. to 300°C. and at pressures up to 1000 Ib./sq. in. 
abs. from the PVT data. Values for the saturated 
liquid, saturated vapor, and superheated vapor are 
presented on a pressure-enthalpy diagram with spe- 
cific volume and entropy as parameters. The entropy 
of the saturated vapor increases with temperature 
over a considerable region, indicating a_ positive 
specific heat of the saturated vapor in contrast to 
the lower alcohols of the same series. 


i em growing importance of processes involving high 
pressure and high temperature conversions has made 
the accurate determination of the pressure-volume- 
temperature (PVT) relationships of many compounds 
increasingly important. As part of an investigation into 
the thermodynamic properties of alcohols, the determina- 
tion of PVT values on normal butanol has been carried 
out to provide data not available previously, and to allow 
compilation and calculation of the enthalpy values in the 
saturated, two-phase, and superheated regions. 


Analytical representation of PVT data is achieved 
through equations of state, of which only a few such as 
those of Benedict, Webb and Rubin “) and of Beattie 
and Bridgman ‘?) have attained practical importance. 
These are all empirical and while theoretical equations of 
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state based on statistical thermodynamics have been at- 
tempted, there has been as yet no dev elopment possible 
of wide application. In this w ork the Beattie and Bridgman 
equation of state has been selected to represent the PVT 
data, and to be used in the calculation of other thermo- 
dynamic properties. 


Graphical presentation of the PVT data and the cal- 
culated thermodynamic functions are also given. 


Available data 


The amount of information on the PVT properties of 
n-butanol is relatively sparse. Vapor pressures up to the 
normal boiling point have been determined by Kahlbaum “) 
and by Butler and co-workers ‘*), Dreisbach (5) has 
reported boiling points at 4 different pressures below 
atmospheric. The definitive normal boiling point deter- 
mination for #7-butanol appears to be that of Wojciechow- 
ski ‘") at 117.726°C., along with a coefficient of increase 
qP _ 9.9372°C. 
dt 
per mm. Hg. Other values accepted by Timmermans (7) 
range from 117.77 to 118.0°C. Tukuzo and Koe ‘*) have 
reported the density of the saturated liquid from 0°C. to 
—80°C. while numerous determinations at single temper- 
atures have been made (7), Critical temperatures have 
been determined several times without reference to other 
critical properties (% 119), and the critical pressure has 
been reported by Herz and Neukirch (2), New values 
have been given recently by Kay and Donham 3) who 
also reported orthobaric densities and vapor pressures up 
to the critical point. The most extensive determination of 
critical properties, including critical density and an equa- 
tion for rectilinear diameters, formed the preliminary 
work in this inv estigation (14), Densities for the saturated 

vapor and the corresponding equilibrium liquids were 
reported at the same time for a temperature range from 
155° to 286°C. A tabulation of critical constants appearing 
earlier (14) has been revised and is shown in Table 1. 


in boiling point with increase in pressure, 


Early work on the heat capacity of m-butanol (1, 1°) 
was limited to providing single values for the liquid, or 
a mean value over a limited temperature range. Parks (17) 
provided low temperature data including liquid butanol 
from its melting point to 20°C., while Williams and 
Daniels ‘'5) gave a value for constant pressure heat capa- 
city as a function of temperature for the range 30-80°C. 
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TABLE 1 
CRITICAL PROPERTIES OF NORMAL BUTANOL 








Date} t°C | P. atm. Pc g./C.C. | Investigators 

1883 | 287 4 | | Pawlowski ay 

1884 | 270.5 | de Heen ‘'® 

1923 | | 48.4 | Herz and Neukirch 

1943 | 288.0 | | Fischer and Reichel ‘? 

1955 89.74 43.6 | 0.267 | Kay and Donham “'® 
48.60 | 0.2700 | Singh and Shemilt “” 


sa 
UW 
| do 
| oo 
a 
oO 
Vi 


Francis 1%) assembled previous data in the form of tem- 
perature dependent equations. 

The only values for heat capacities in the vapor state 
were given by Sinke and de Vries ‘2°) for the range 120- 
165°C. at 750 mm. pressure. 


Experimental 
1. Materials 


The n-butanol (c.p. grade) was distilled at high reflux 
in a glass-helici packed column and middle cuts further 
refluxed and distilled successively from calcium oxide, 
fresh magnesium and activated alumina. In all distillations 
only the middle fraction (one-third) was collected. The 
final product gave a refractive index of 1.3973 at 25°C. 
measured on a carefully calibrated Abbe-type refracto- 
meter, a density of 0. 80565 at 25°C. and a ‘boiling point 
of 117.77°C. using a thermometer calibrated against a 
laboratory standard with an N.B.S. certificate. Two 
separate lots of 7-butanol were prepared and gave identical 
results within the precision of the measurements (0.0001 
refractive index units, 0.00003 in density, and 0.03°C. 
respectively ). 

The mercury used in contact with the 2-butanol was 
triple distilled. ‘Chlorobenzene, aniline, methyl salicylate, 
bromonaphthalene, n. butyl phthalate, and benzophenone 
used for the constant temperature column were reagent 
grade or better and were further purified by distillation 
at high reflux. 


2. Apparatus and procedures 


The apparatus, designed after that of Kay (!, 22) con- 
sisted of an adjustable steel block filled with mercury and 
connected to a calibrated glass capillary containing the 
sample and enclosed in a constant temperature apparatus. 
The pressure in the capillary of 2 mm. true bore tubing, 
was transmitted by mercury to an oil line leading to a 
dead weight gauge tester. ‘A piston in the compressor 
block could be moved to adjust the pressure and keep 
the mercury-oil level constant, as shown by an electrical 
relay arrangement. For temperature measurement a 
oo -constantan thermocouple was used with a type 

2 potentiometer having a Weston Standard cell w ith an 
N. BS. certificate. The temperature was kept constant 
around the capillary tube by boiling a purified organic 
compound under a controlled pressure. By selection of 
the proper compound a wide range of temperatures was 
obtainable. The pressure was measured by a manometer 
and controlled by a manostat operating through an elec- 
tronic relay. 


The constant temperature vapor bath was a silvered, 
vacuum-jacketed distillation column 4 feet long with 
sufficient expansion bellows to allow a temperature dif- 
ferential up to 340°C. 
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The compressor block contained the mercury-oil 
interface detector. The capillary was mounted by packing 
and a special-type gland. The piston was sealed with a 
special V- leather packing and could be moved in or out 
without loss of mercury up to operating pressures of 
1200 lbs./sq. in 


The internal diameter of the “constant-bore” capillary 
tubing throughout its total length, and the volumes in 
the region of tubing joins and the volume of the sealed 
tip were obtained by weighing amounts of triple distilled 
mercury released from the tubing and measuring the 
height difference with a Gaertner cathetometer previously 

calibrated against a metre stick with an N.B.S. certificate. 


The capillary tubing was filled by means of a carefully 
cleaned filling tube to a depth of about 6 inches with 
purified #-butanol which had just been boiled vigorously 
to remove dissolved gases. The tubing was inserted into a 
simple filling apparatus where a vacuum (3 x 10° cms. of 
mercury pressure, was applied until about 50% of the 
n-butanol was vaporized and removed. Purified mercury 
contained in the filling apparatus was then spilled into the 
capillary tubing sealing in the liquid butanol sample, and 
the tubing was inserted into position in the pressure block. 


The copper-constantan thermocouple was constructed 
from Leeds and Northrup glass-insulated duplex thermo- 
couple wire No. 30 B. & S. gauge, and was calibrated 
against a 2.5 ohm platinum-resistance thermometer with 
an N.B.S. “coteilicete: The latter was also independently 
checked against the ice-point. The av erage deviation 
throughout the range to be used was + 0.03°C. The dead 
weight test gauge, manufactured by Barrett, had a range 
up to 4000 p.s.i. Piston diameters and the weights used 
were checked for accuracy and found to agree with the 
manufacturer’s listing (i.e. maximum total error less than 
0.1%). 


The general operating procedure involved manipula- 
tion of the compressor block piston and the weights on 
the test gauge pan to maintain correct mercury levels. 
Beginning with a liquid sample, vaporization was carried 
out to give as large a vapor-to-liquid ratio as possible for 
all two-phase readings. Condensation from a vapor state 
was also carried out to ensure the approach to equilibrium 
from “both sides”. Thermal equilibrium was considered 
established when temperature readings changed less than 
0.1°C. for a period of one hour, and the mercury level 
changed less than 0.5 cm. in the same period, indicating 
paces no vaporization or condensation and hence no 
large heat transfer demand. For this period also, room 
temperatures were maintained constant to 1°C. and 
atmospheric pressure continually noted. Following estab- 
lishment of equilibrium, readings were made of the 
weights on the gauge tester, and of mercury levels 
throughout the system including meniscus height, and 
position of liquid meniscus when present. Temperatures 
ranged from 117.6°C. to 305.5°C. Below the critical point 
both liquid and superheated vapor systems as well as the 
two-phase equilibrium system were examined. 

The actual weight of the samples used was obtained 
by measurements of liquid volume below 80°C. and utiliza- 
tion of known liquid density data ‘7). Isothermal volume 
measurements for a completely and partially vaporized 
sample were used to give saturated liquid and vapor 
densities. Dew points were obtained by increasing pressure 
incrementally on the superheated vapor at constant tem- 
perature, and noting the point when the superheated 
vapor specific volume decreased significantly. 

Two completely independent samples were used for 
purification and for measurement purposes. Comparison 
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of boiling-points, refractive indices, and vapor pressures 
at two different temperatures showed no difference 
beyond that inherent in the methods of measurement. 

Corrections were applied for the mercury meniscus in 
accordance with values given in the International Critical 
Tables (23). Volume corrections for the expansion of the 
glass capillary tubing due to temperature were made 
according to the equation ‘*4), 


V = V,(1 + at) 
where a = cubical expansion coefficient, 9.9 x 10°, 
and t = temperature change in °C. 


The effect of pressure on the glass capillary volume was 
calculated Y the equation 25) 
AV ’ (30.27 x 10 'p + 99x 10-12?) 
“ pa ote p = pressure, Kg/cm? 
and V,, = original volume, cc’s. 


Results and calculations 
1. PVT Measurements. 

The results of the PVT measurements are presented in 
Table 2 and graphically in Figure 1, and cover a range of 
temperatures from 197°C. to 305°C. and pressures from 
100 to 1000 psia. 

The original results obtained as isothermal data were 
transformed to isometrics by a method similar to that 
used by Andrews ‘"), and modified by Beattie and Bridg- 
man (2), The transformation was effected by calculating 
the quantity oe -RT) for each measurement on a given 
isothermal. A straight line was then passed through the 
point of zero density (where PV-RT = O), and the 
highest measured density on the isotherm. Deviation of 
the observed values from this line were plotted against 
density. From this plot values of PV-RT and hence of 
the pressure (considered as observed pressure) were inter- 
polated for a set of evenly spaced densities, the same 
density being used at each isothermal. The data, smoothed 
with respect to volume are presented in Table 3. 


Vapor pressure results ranging from the normal boil- 
ing point to the critical point are given in Table 4, along 


with saturated vapor and liquid densities previously 
reported 14), 
2. Equation of state. 
The PVT data were used for calculation of the 
constants in the Beattie-Bridgman equation of state () 
where 
RVC k \ 
P = (V +B 
V2 \ 


V1 


The isometrics (Table 3) were each straightened by using 
factor C, and that value of C was chosen 
which straightened out the most curved isometric. The 
Beattic Bridgman procedure ‘*) was followed using a 
‘corrected pressure” p} defined as 
p' = p+ F(V,1 ryv) — dV) 

For each isometric, plotted as p! vs. T, there were ob- 
tained values of W(V) and @(V) appearing as slope and 


a correction 


intercept Since 
\ 1B V2, W\ ) 
Kk 
and A V2, d(\ ) 
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TABLE 2 


ISOTHERMAL PRESSURE-VOLUME DATA (OBSERVED) 
FOR n-BUTANOL 





P = Pressure, p.s.i.a. V = Volume, ccs/gm. 
Temp = 197.28°C. | Temp. = 268.90°C. | Temp. = 286.99°C, 
P V P V P | Vv 
£13.72 59.63 138.90 60.43 405.31 17.40 
122.87 54.28 457.65 | $2.53 454.63 | 14.60 
127.49 51.84 195.84 40.92 503.92 | 12.00 


128 


“HA =~ 00 | 


94 | 49.72 244.33 | 31.34 602.74 | 10.20 


293.29 | 24.86 632.41 | 7.10 
is _ | 342.53 | 20.00 661.86 | 5.85 
remp. = 209.11°C. | 393.93 | 16.10 | 671.60 | 5.50 
<canticiedeiic ise 13.00 711.28 3.170 
; 481.03 10.70 731.06 2.590 
P \ 751.03 | 2.403 
ae | 800.97 2.300 
73 | 59.83 850.96 | 2.223 
26 | 50.36 | Temp. = 272.25°C. | 900.96 | 2.203 
07 42.80 cat 1000.94 | 2.195 
19 37.87 


ee ae 
Sow 
~ 


Temp. = 221.80°C. 


Pp | v 


278.13 28.25 


a 336.38 21.30 | Temp. = 294.68°C. 








p Vv 355.87 19.40 eS 
404.87 | 15.70 
: ‘ 
mn | ee |) See | ee 
oc, | coon | eee | Ole. | Set | ome 
180.27 37 69 650.85 | 2.218 | 405.68 | 17.56 
99 51 37.95 | 700.84 2.170 | 454.94 | 15.53 
214 02 29 84 800. 82 2.110 504.19 | 13.44 
900. 80 2.051 553.36 | 11.52 
1000.79 2.014 603.09 | 9.76 
Temp. = 232.45°C. 652.56 | 8.13 
E 702.28 | 6.73 
| 751.17 5.674 
P \ Temp. = 276.48°C, | 801.22 4.863 
| 851.06 | 4.296 
124.06 61.99 | 900.99 3.850 
142.62 52.43 P Vv | 950.96 3.459 
161.50 45.34 | 1000.45 | 3.053 
180.63 39.54 278.04 28.37 
199.83 34.75 355.79 119.59 
hi ree m4 a 404.87 | 16.14 
£99.09 4.94 453.96 | 13.37 | Temp. = 300.56°C. 
258.11 24.31 552.57 8.646 oi 
263 91 | 23.34 601.21 | 7.160 | hon. 
204.7 23.32 | 700.55 2.180 P V 
800.52 2.117 - | ——- 
: o¢ | 900.50 | 2.083 | 308.09 | 27.29 
Temp. = 244.63°C. | 1999" 48 2.040 | 366.56 | 22.36 
405.77 | 18.31 
p V 454.97 | 15.67 
Temp. = 283.94°C. 504.30 13.64 
129.01 61.35 |_ ati a 2s 3 
147.64 52.64 . co ‘= 
166.58 | 45.76 Pe 652.79 | 8.64 
oo | ae 702.41 | 7.29 
204.98 | 35.61 
243.83 | 28.38 | 288.16 | 28.37 | 752.01 | 6.28 
282.98 23 02 356.25 | 20.52 801.66 | 5.472 
302 61 20.69 405.23 | 17.30 851.31 | 4.782 
ae : 901.09 4.255 
322.29 18.63 | 454.40 | 14.30 ; 7 
327.24 18.13. | 503.98 | 12.10 | 950.90 3. 850 
553.01 9 92 1000.02 3.513 
602.49 7.83 | 
Temp. = 260.55°C. | 651.48 | 5.80 | 
730.86 | 2.439 | Temp. = 305.52°C. 
800.82 | 2.297 | 
P V 900.79 | 2.220 
1000.77 | = 2.156 | , V 
134.05 | 61.44 | 
152.74 53.09 279.26 | 31.07 
171.52 | 46.45 | 308.23 | 27.68 
210.16 36.59 | 376.32 | 21.04 
258.83 28.20 | 405.92 | 19.68 
307.88 22.19 455.10 | 16.89 
357.17 17.73 | 504.43 14.61 
406.57 13.94 | 603. 86 12.68 
421.41 12.93 653.34 9.68 
703.34 8.34 
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Figure 1—Pressure-volume isotherms for n-butanol. 


the values of B and A for each isometric were found, and 
plotted as functions of density. These plots enabled the 
four remaining constants in the equation of state to be 
determined and are shown in Figures 2 and 3. A definite 
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Figure 2—The Beattie-Bridgman function, A = V2  (V). 
Related to density for n-butanol. 


break occurred in the two plots at the region of p = 3.00; 
and hence two different sets of constants were obtained. 
Table 5 presents these compared with those of methanol, 
the only other normal alcohol for which results are avail- 
able (27), 
By using the virial form of the equation of state 
P=RTp+6.p?+y7.p? +6. pt, 

the virial coefficients may be evaluated from the Beattie- 
Bridgman constants as, for example, with the important 
second viria coefficient, 


§ = B, — — = 


3. Compressibility 
T 


P 

RT’ 
from the experimental PVT data for isotherms from 
197°C. to 305°C. and the values presented in Table 6. 
These results, expressed as a function of reduced pressure, 
are also given in a “compressibility chart’, Figure 4, with 


The compressibility factor, Z = was evaluated 


TABLE 3 


ISOMETRIC PRESSURE-TEMPERATURE DATA (OBSERVED) FOR n-BUTANOI 


Absolute pressure 


Density 
gm. mole/(p) | 

litre 260 265 270 27 

“0.50 18.00 18.42 | 18.85 | 19 
0.75 23.60 24.11 24.70 25. 
1.00 27.90 28.70 29.60 30. 
1.25 30.40 31.50 32.65 | 33. 
1.50 32.50 33.90 | 35.30 | 36. 
1.75 33.35 35.00 | 36.70 | 38. 
2.00 | | | | 
2.50 | 
3.00 | 
3.50 | 
4.00 | 
4.50 
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80 
40 


atmos. 
Temp. °C. 


| 
280 | 


285 200 295 300 
19.60 20.00 20.40 20.80 21.20 | 
25.00 | 26.37 26.95 27.50 28.05 
31.23 | 32.10 32.95 33.80 34.65 
35.00 | 36.15 37.30 38.50 39.68 
38.25 | 39.70 41.20 42.65 44.10 
|} 40.13 | 41.82 13.58 45.30 47.00 
| 41.50 | 43.40 45.45 47.50 49.50 
44.30 | 46.90 49.55 52.20 54.80 
47.60 50.50 53.40 56.40 59.20 
49.80 53.50 57.20 60.95 64.60 
52.50 | 56.70 60.85 65.00 69.20 
54.70 | 59.58 64.40 69.20 74.00 
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PABLE 4 
POR PRESSURES AND ORTHOBARIC DENSITIES FOR n-BUTANOL 
density G/CC 

€ ( Pre ‘ = 

Liquid Vapor 
117.58 4.65 
12¢ RS 20 36 
131.835 23.79 
155.77 48 61 0.70520 0.0079 
155.99 48.74 
167.74 66.97 
173.30 77.17 9 67059 0.01249 
175.78 81.64 
183 ¢ 97 8&5 0.62478 0.01587 
198.89 137.17 ). 61920 0.02106 
209 11 170 19 
212.02 179 16 0. 59305 0.02775 
221.59 215.8) 
221.80 216 53 
232.45 264 41 0.04280 
237 .31 287 .06 0 53590 
244.63 326.99 Q 51993 0.05516 
260 55 421.22 () 48589 0.07733 
268 90 4&1 03 0.45740 (09480 
272.2 541 .0* 0. 43200 0.1023 
276.5 601 0* 0.42370 0.1300 
283.9 653 .0* 0.3620 0.1681 
285.7 0 .3220* 0.2402 

ext polated values 


temperature as a parameter, and the calculated values in 


the form of smoothed curves. 


4. Volume, vapor pressure, and heat of vaporization. 


Saturated lig ud volumes up to 100°C, were calculated 


from the Bridgman equation ‘*°) 


dG 0.82390 0.699 10-4 + 0.32 K 10°%2, (1) 


where t temperature oy A 

These values on a Jarge-scale plot, with the experimental 
| 

ones from 155°C. to 287 ¢ illowed a smooth interpola- 

tion for the range 100° to 155 C. Saturated vapor volumes 

up to 155°C. were calculated using the compressibility 


factor from the generalized charts of Hougen and Watson 


The Vapol pressure values available were plotted on 


a large scale drawn 
through the 


c 


Cox-type graph, a smooth curve 


data, and apor pressures read for the 


desired increments of temperature 
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Figure 3—The Beattie-Bridgman function, 
B = V? ¢ (V)—V. 


Related to density for n-butanol. 


Latent heats of vaporization were calculated by the 
well-known Clapeyron equation (), by the Othmer 
method ‘), and by Othmer’s method as modified by 
‘Thomas (42), The values obtained by the three methods 
varied by less than 3% in all cases. A smooth curve 
averaging the results and extrapolating to a value of zero 
at the critical temperature was used to portray the best 
value. 


5. Enthalpy 

Liquid and vapor enthalpies were obtained from the 
enthalpies of the ideal gas, the equation of state, and the 
heats of vaporization in the following manner. First, the 
enthalpy of the ideal gas at zero pressure was obtained 
from a heat capacity equation. For this the heat capacity 
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IS FOR BEATTIE-BRIDGMAN EQUATION OF STATI 


n. butanol methanol 


400 250) 300 
4.00 4.00-4.50 
0.10730 0.09246 


14.245 $43. 409 
0.1047 0.09929 
0.6571 0.60362 

10% 1.00 K 10° $.2031 &K 10° 
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Figure 4—Compressibility factors for n-butanol. Figure 5——Pressure-enthalpy diagram for n-butanol. 
TABLE 6 (a) 
COMPRESSIBILITY FACTORS FOR n-BUTANOL 
Volume — litres/gm. mole Press. — atm. Temp. — °C. 
es 5‘ ‘i ee - a hate ‘a a j 
Temp. | Volume 1.50 1.75 2.00 2:25 Z00 | 2.85 3.00 3.25 3.50 4.00 4.50 
197.3 | Press. | | 10.23 9.79 9.27 8.39 | 7.65 
Z | | 0.7944 0.8240 0.8403 0.8692 0.8915 
209.1 | Press. 12.28 | 11.63 11.00 10.38 9.78 8.77 7.95 
e 0.7758 0.8082 0.8339 0.8525 0.8650 0.8865 0.9090 
221.8 | Press 15.58 14.38 13.33 12.42 11.62 10.88 10.23 9.16 8.23 
Z 0.7684 0.7979 0.8220 0.8423 0.8596 0.8720 0.8838 0.9041 0.9132 
232.4 Press 17.78 16.27 15.00 13.87 12.90 12.04 10.227 10.60 9.950 8.580 
: 7 0.7499 0.7845 0.8133 0.8357 0.8549 0.8726 0.8835 0.8941 0.9156 0.9305 
y 
S 244.6 | Press. 20.75 18.92 17.20 15.71 14.50 13.45 12.56 11.74 11.03 9.85 8.86 
> z 0.7261 0.7726 0.8025 0.8249 0.8457 0.8632 0.8794 0.8901 0.9006 0.9192 0.9290] 
) — —— — 
; 
ABLE 6 (b) 
, COMPRESSIBILITY FACTORS FOR n-BUTANOI 
, Volume litres/gm. mole Press. — atm. emp Cc 
, ; | 
Temp. Volume, 0.222 0.25 0.2857 0.333 0.40 0.50 | 0 5714 0.666 0.80 1.00 1.33 2.00 
260 | Press. | 133.35 32.50 | 30.40 27.90 23.60 18.00 
Z | | 0.4356} 0.4952 | 0.5559 0.6377 0.7190 0.8228 
| | 
265 | Press 35.00 33.90 31.50 28.70 24.11 18.42 
; 0.4526) 0.5109 0.5703 0.6495 0.7273 0.8339 
| 
270 | Press | 36.70 | 35.30 32.05 29.00) 24.70 18.85 
’ | 0.4705) 0.5275 0.5800 0.0641 0.7387 0.8459 
1 
| 
275 Press 38.40 30.80 $4.85 40.40 25.25 19>) 
Z | 0.4878) 0.5449 0.6020 0.6758 0.7483 0.8548 
| 
280 Press 54.70 S52. 50 19.80 | 47.60 44.30 41.50 10.18 38.25 35.00 | 31.23 "5 80 19 60 
Z 0.2675 0.2891) 0.3134) 0.3492) 0.3904) 0.4571 0.5015} O.5617 | O.6108 | O. 6880 0.7516 0.8036 
| | 
285 | Press. | 59.58 | 56.70 | 53.50 | 50.50 | 46.90 | 43.40 | 41.82 | 39.70 | 36.15 32.10 6.3 0.00 
Z 0.2890) 0.3095) 0.3377) 0.3675) 0.4096) 0.4738) 0.5218) 0.5778 0.0314 0.7008 0.7675 0.8736 
290 Press 64.40 60.85 57.20 53.40 49.55 45.45 43.58 41.20 37.30 32 95 26.05 20.40 
Z 0.3096; 0.3291 0.3536) 0.3851 0.4289) 0.4918) O.5388) 0.5943 0.0457 0.7130 0 } 0.8829 
295 Press 69.20 65.00 60.95 56.40 52.20 47.50 $5.30 42.05 38.50 33.80 7 50 0 80 
Z 0.3298 0.3485) 0.3735) 0.4032) 0.4478 0.5094) 0.5552) 0.6093 0.0000 0.7250 0.8020 0.8923 
? 300 Press 74.00 69.20 64.60 59 20 54.80 49 50 47.00 14.10 349 OS 4.05 Ss US 21.20 
Z 0.3493 0.3678 0.3924 0.4192, 0.4660; 0.5262) 0.5710; 0.6245 0.6749 0.7360 0.7950 0.9015 
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of the vapor at zero pressure was determined by using the 
characteristic frequencies of all the intramolecular bonds 


as measured spectroscopic: ally, 


(Cp)o = — 8.324 X 10-3 


where T = temperature in °K. 


and calculating by means 
of the appropriate Einstein functions. The frequency 
values of Bennewitz and Rossner ‘%), as subsequently 
modified by Dobratz 4), were used. The heat capacity 
equation thus obtained was 


> + 1.2778 X 10-*T — 0.48048 X 10-°T?, 


(2) 


Changes in enthalpy with pressure and temperature may 
be expressed by the equation 


dH = 


CAT + k ey (; 


ON 


oT 


> 


\Je. 


(3) 


By integration of this equation and selection of an appro- 
priate standard state, enthalpy values for the non-ideal 
saturated vapor state were obtained. The standard state 
selected was a value of zero enthalpy for the liquid at 
0°C. The temperature term was integrated using Equation 
The pressure term was integrated by the 
method of Hougen and Watson ‘%5) up to 250 °C. For 


(2) above. 


the range 250 


were 


to 300°C., 
determined from 


the 


constants 


values of the pressure integral 
of the Beattie- 


Bridgman equation of state, as suggested by Dodge (36) 
using a differentiation of the expanded form “of the equa- 
tion of state. The enthalpy value for the ideal gas at 0°C. 
and zero pressure was calculated by adding the heat of 
vaporization value at 0°C. to the standard state value for 


the liquid. 


Liquid enthalpy values were then calculated by addi- 
tion of latent heats of vaporization to the enthalpy values 
y were also 
determined by Williams and Daniels’ equation (1%) 


for the saturated vapor. 


G 


= 0.49992 


where t = temperature in °C. 
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THERMODYNAMIC PROPERTIES OF n-BUTANOL 


Volume 
liquid 


CCs. 


pd OD IND ID meet net ee ee et et ee et eet eet eek ee ee et pee pee eet tee fee fee fee pee fe fee fee fet 
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.2137 
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ey 


+ 0.00263t + 0.0000018t2, 


me 
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(4) 


Enthalpy 


Up to 70°C. the values checked within 0.5% 
6. Entropy. 


Entropies were calculated in a manner similar to that 
for enthalpies. Values for the vapor were obtained from 
the equation (37) 


where P, 


0, 


T 


a s= s+ [o 


oO 


P 
ver] 
-T 
Po 
ees 


O\ 
oT 


) SHEE ab eee (5) 


and S, = an arbitrary standard state value of entropy 
at 0°C. and zero pressure, and taken to be 
zero for liquid. 


The temperature integral was evaluated using Equation 
The pressure integral was determined for the range 
of 250°C.-300°C. from the Beattie-Bridgman constants as 
before, and up to 250° by the method of Hougen and 
The values at 250°C. 
checked within 1.5%. 


(2) 


Watson (38), 


by the two methods 


Entropies for the liquid were obtained from the 


equation 


where S,, 


and T = temperature in °K. 


Si 


- A 


= - 


H, 


T 


. (6) 


entropy of liquid, cal. per mole per °K. 
entropy of vapor, cal per mole per °K. 
latent heat vaporization, cal. per mole. 


Liquid entropies were also determined for the temper- 
ature range 0°C.-80°C. from Daniel’s equation of heat 


capacity and were found to agree within 0.8%. 


Enthalpies and entropies for the superheated region 
were determined in a similar fashion to the saturated 
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SATURATED ZONE 
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vapor by using the Beattie-Bridgman constants from 
250°C.-300°C., and the Hougen and Watson method (89) 
from 117.7°C.-250°C. Values below 117.7°C. were obtained 
by extrapolation. 


- 


7. Thermodynamic network. 


The thermodynamic properties measured and calcu- 
lated as above are presented for the saturated zone in 
Table 7, for the superheated region in Table 8, and 
graphically in Figure 5. 


Discussion of results 


A high degree of purity for the n-butanol was ob- 
tained. The density measurement of 0.80565 compares to 
the values of 0.80567, 0.80572, 0.8057, 0.8057 selected by 
Timmermans (7), The refractive index value of 1.3973 
compares to the critically selected values of 1.3974 and 
1.3970. The normal boiling point obtained during purifi- 
cation was 117.77 + 0.03°C. After the purified material 
was enclosed in the PVT apparatus, a series of readings 
at different times gave an average value of 117.69 + 0.03°C. 
This agreement between a static and dynamic method 
appears quite acceptable. It also compares favorably to 
Timmermans (7) selection of 10 values ranging from 117.7 
to 118.0°C., and Kay and Donham ‘!%) who reported 
117.73°C. 


In addition a further criterion of purity was used. 
The average pressure rise on condensation was deter- 
mined from the dew-point pressure and the pressure 
value at complete condensation. Over a series of temper- 
atures from 220°C. to 270°C. this value averaged 0.15 p.s.i., 
somewhat more favorable than the single value of 0.32 
p.s.i. reported by Kay and Donham “%) at 230°C. A 
similar gauge of the purity showed in determinations of 
vapor pressures with widely varying percentages of 
vaporized material. At 185° C.. for example, the deviation 
was less than 0.04% for a change from 20% to 80% 
vaporized. 

All the usual corrections to pressure and volume were 
applied, with the exception of any correction for the 


presence of mercury vapor. The usual correction consists 
of subtracting the vapor pressure of mercury from the 
measured pressure, allowing for the effect of the total 
pressure on this vapor pressure value (the so-called 
Poynting effect). Jepson and Rowlinson ‘*) have recently 
discussed this correction and applied additional corrections 
for deviations from ideality in the gaseous solution of 
mercury and the contained vapor. However, in the 
measurements reported here an extremely small mercury 
surface was available to the system. In addition precise 
duplication of results on the two- -phase system when 
approaching it both from a superheated vapor and a total 
liquid state, confirmed that no mercury vapor correction 
was required. 


The vapor pressure data reported here show complete 
consistency with Butler’s data below the normal boiling 
point when plotted on a large scale Cox-type chart. 
Moreover both an Antoine equation and another type of 
a 3-constant vapor pressure equation show no discon- 
tinuities at the boiling point when applied to all available 
data (41), The isothermal plot (Figure 1) and the calcu- 
lations used in determining the critical points 4) both 
indicate a high degree of ‘self- consistency of the results 
given here. It is thus considered that precision of the 
measurements (0.01 atmospheres for pressure; 0.03°C. for 
temperature; and 0.5% for volumes) is satisfactorily re- 
flected in the PVT values reported. 


While data in the superheated region are available 
only from 197°C. to 305°C., measurements on properties 
of the saturated vapor and liquid are given from the 
normal boiling point to 287°C. The size of the sample 
that could be put in the apparatus limited the range of 
vapor volumes that could be measured at the lower 
temperatures. The fact that measurements could be made 
above the critical point meant however that the significant 
range of values was adequately covered. 


The vapor pressure data used for the presentation of 
the thermodynamic network are based on actual measure- 
ments throughout the temperature range from 0°C. to the 


TABLE 8 


THERMODYNAMIC PROPERTIES 


OF n-BUTANOL 


SUPERHEATED REGION 


Temp. — °C. 











Press. — p.s.i.a. 


volume, V — cc's/gm. = 62.43 (cu.ft./Ib.) 
1 esi nes 
enthalpy, H — cal./gm. = (B.t.u./Ib.) (; .) Entropy, S — cal./(gm.) (°C.) = B.t.u./(Ib.) (°F.). 
Temp. 
Press. pe ; . ¥ “a | q 
117.6 120 130 140 150 160 170 180 190 200 
14.7 V\ 428.37 425.26 440.00 451.94 463.11 | 474.54 485.49 496.70 507 .92 | 519.15 
H| 213.89 214.88 219.20 223.48 227.97 232.58 237.16 241.87 246.64 251.50 
S 0.5596 0.5629 0.5780 0.5847 0.5955 0.6061 0.6167 0.6274 0.6378 0.6481 
25.0 \ 264.07 270.73 277.01 | 283.79 289.69 297 .02 304.79 
H 222.85 227.48 232.19 | 236.57 241.25 246.07 250.98 
S 0.5691 0.5799 | 0.5907 | 0.6012} 0.6122 0.6241 | 0.6331 
| | 
50.0 \ 134.33 | 137.27 140.68 143.78 148.13 
H 230.23 | 234.84 239.60 244.40 249.51 
S 0.5691 0.5802 0.5911 0.6020 0.6128 
750 V 91.82 95.14 97.21 
H 238.01 243.04 247.95 
S 0.5776 0.4889 0.5999 
100 V 69.63 tag 
H 241.39 | 246.32 
S | 0.5785 0.5898 
| 
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critical point. The results of Kahlbaum ‘) and of Butler 
(4) up to the normal boiling point agree very well. From 
there to the critical point the results reported here are 
considerably higher than those reported by Kay and 
Donham (13), The difference is approximately 2% until 
close to the upper limit when it increases to 10%. There 
appears to be no obvious reason for the disagreement. 
It is felt, however, since a significantly lower pressure 
rise on condensation was obtained, and extreme caution 
used in establishing of equilibrium, that the values reported 
here are more reliable. One other interesting fact is that 
the application to the normal alcohols of the Riedel ‘*) 
and Lydersen ‘43)group contribution method for estima- 
ting critical properties, with due allowance for increasing 
chain length and molecular weight, indicates that the 
results reported here are more consistent with those of the 
other normal alcohols (41), Better agreement with an 
Othmer-type of reference substance relationship is also 
obtained, and in addition a large scale plot of log p against 
1/T shows the slight tendency to an S-shaped curve that 
has been noted for the corresponding paraffins (#4), These 
results were therefore used for all the thermodynamic 
calculations given here. Stull ‘*5) has previously calculated 
vapor pressures from the normal boiling point to the 
critical point by an unstated extrapolation method. These 
deviate from the measured values used here by an average 
of 2%. 

Evaluation of the constants for the Beattie-Bridgman 
equation of state was limited to the temperature range 
of 250-300°C. because insufficient volumetric data were 
available at lower temperatures. While such an equation, 
with 5 empirical constants, has been shown to fit PVT 
relations with a high degree of accuracy, its adequacy is 
illustrated by the necessity of using two sets of constants 
to cover the range of densities in this case. The values of 
the constants obtained for 7-butanol are similar in magni- 
tude to those for methanol and indicate their possible 
usefulness in homologous series calculations. They have 
also been applied to calculate other thermodynamic pro- 
perties, and are of interest in calculating virial coefficients 
and hence for comparison with theoretical calculations of 
intermolecular forces and force constants. 


The specific volumes of the superheated vapor at 
temperatures below 200°C. were not of as high order of 
accuracy due to the size of sample that had to be used, 
and the ‘difficulty of actual measurements of large volumes 
at low temperatures and pressures. As a result the use of 
a generalized method appears satisfactory at temperatures 
under 200°C. At higher temperatures and particularly at 
higher pressures significant deviations occur from the 
usual generalized compressibility chart. From the PVT 
measurements reported here, the calculation of compres- 
sibilty factors has resulted in a chart (Figure 4) of limited 
but useful range. Some of the compressibility factors vary 
as much as 40% from those predicted by the gener ralized 
compressibility chart. Below the critical point the correc- 
tions are almost all negativ e, and above the critical point 
are all positive. The increasing use of a third parameter, 
such as Pitzer’s acentric factor | 46) or the critical com- 
pressibility (47), enables much more accurate generalized 
predictions. For intensive application such as in an homo- 
logous series, this still does not attain the accuracy possible 
through use of experimental data. 


The close agreement shown between liquid enthalpy 
values based on Daniel’s equation (measured heat capacity 
values ) and those obtained by calculation from spectro- 
scopic data and addition of heats of vaporization, indicates 
a high degree of consistency. It also lends validity to the 
use of spectroscopic calculations at higher temperatures 
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where liquid heat capacity values were not available. 
While the vapor heat capacity values of Sinke and de 
Vries (20) are approximately 10% higher than those cal- 
culated with the data of Bennewitz and Rossner 33), the 
former were actually limited to a determination af 5 
values over a narrow range. The overall effect on enthalpy 
values would be less than 0.3% in any case, and hence 
the use of the full range of calculated values was con- 
sidered to be justified. Since the entropy values as calcu- 
lated were based on the same data, the good accuracy of 
these is also indicated. 


The thermodynamic properties reported are in reality 
overall or apparent values and do not take into account 
the effect of association between 7-butanol molecules. 
From the practical standpoint such apparent values repre- 
sent the actual behavior of n-butanol independent of 
whether or not the state of equilibrium between associated 
and unassociated molecules is known. 


The entropy values for the saturated vapor are inter- 
esting from the design standpoint since there is a con- 
siderable region through which positive specific heats are 
indicated. This means that adiabatic compression from 
certain superheated vapor regions leads to the compressed 
vapor becoming “wet”. In contrast to this both methanol 
(8) and ethanol (**) have saturated vapor entropies that 
decrease continuously with increasing temperature. The 
normal paraffins exhibit similar behavior (5°), 


From overall considerations of accuracy, it is estimated 
that the error in the thermodynamic properties now pre- 
sented for -butanol is within 1% for most of the range 
covered and not much in excess of that even at highest 
temperatures and pressures. 


Nomenclature 

a, A, Ao, b, B, Bo, C, E = constants in Beattie-Bridgman equation 
of state. 

C, = heat capacity, constant pressure 

d = density (Bridgman equation) 

d = differential operator 

H = enthalpy 

AH, = latent heat of vaporization 

p! = corrected pressure in Beattie-Bridgman equations 

p = pressure, kg./cm.? 

r = pressure, psia. 

R = universal gas constant 

S = entropy 

t = temperature, °C. 

a = absolute temperature 

V = specific volume 

AV = volume change 

z = compressibility factor 

Greek letters 

a = cubical expansion coefficient 

B, y, 6 = virial coefficients 


p = density 
A = finite difference 
@ = function 


= function. 


Subscripts 

L = liquid 

O = zero (pressure, temperature ); original value (volume). 
; vapor. 
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Rate Control in the Hydrometallurgical 
Preparation of Uranium Dioxide’ 


D. R. WILES? 


Recent studies of the hydrogen-precipitation of 
uranium dioxide from alkaline carbonate solutions 
have been extended. It was revealed that the reduc- 
tion takes place in two phases. The first of these is 
the one described by Forward and Halpern. The 
second stage, much slower than the first, appears to 
be associated with the complete covering of the 
catalyst by the solid product. It is shown that while 
quantitative information is unreliable, a good quali- 
tative description of the second-stage reaction can 
be based on the diffusion of reactants through this 
product layer to the catalyst surface. 


NE of the more active areas of research in hydro- 
C metallurgy in recent years has been that of the 
oe ee of metals from solution of hydrogen reduc- 
tion (1. 2.8), By extension of this method, Forward and 

Halpern ‘4 found it agen to strip uranium from car- 
bonate leach solutions by hydrogen precipitation of uran- 
ium dioxide. Analogous reactions have been investigated 
in vanadium (5). A solid catalyst, such as nickel powder, 
is used. 

The reduction of uranium was found ‘4) to occur as 
described in Equation (1): 

(Ni) 


UO(CO;)s. + Hy + UO. + 2HCO;- + CO;" (1) 


Forward and Halpern also studied the kinetics of the 
reaction. They found that the rate of the reaction is 
proportional to the catalyst surface area and to the first 
power of the hydrogen pressure, but that, except at the 
very end of the reaction, the rate is independent of the 
uranium concentration. They thus found that the rate of 
the reduction could be described by the equation, 


d{U] 


k{Nij. [H»2] e#/k! ( 
dt 


nN 


During the course of a further investigation of this 
process, it has been discovered that the precipitation 
proceeds in two phases, giving the same product, but 
with quite different rates. It was found that, following 
the stage which was investig: ated by Forw a = Halpern, 
there occurs a slower st: age, char acterized, appears, by 
the fact that the catalyst is completely ati. with a 
layer of the reaction product. It is likely that this con 
dition will obtain in any industrial applic ation of the 
process. 


i\Manuseript received March 15, 1959 
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Bunji and Zogovic have recently shown (6) that the 
product UO, can itself function as a catalyst for the 
reduction. However, catalysis by uranium oxide results 
in reduction at a rate very much slower than even the 
slow second stage in the presence of nickel. 


Experimental 


The experiments to be described were performed using 
solutions of uranium trioxide in ethylenediaminium car- 
bonate. This cation offers two advantages in the present 
case: its solutions are more stable than those containing 
ammonia, and the reduction rates obtainable in the amine 
solutions are much higher and more convenient for study 
than those obtained in sodium carbonate solutions. Since 
the cation does not appear explicitly in E quations (1) and 
(2), one would expect the choice of cation not to have 
a qualitative effect on the reaction. However, it will be 
pointed out later that there is a substantial difference in 
the rates in solutions using different cations. 


Solutions were prepared by passing gaseous carbon 
dioxide into suspensions of uranium trioxide in aqueous 
ethylenediamine. After dissolution of the uranium, the pH 
was adjusted to 9.5 by addition of further CO,. The 
uranium trioxide was prepared by heating uranium per- 
oxide which was, in turn, precipitated with hydrogen 
peroxide from a sulphuric acid solution of sodium 
diuranate. 

The ethylenediamine was technical grade, made by 
Carbide and Carbon Chemicals Corpor ation. Carbon 
dioxide and hydrogen were obtained in cylinders from 
liquid Carbonic Canadian Limited, and Canadian Liquid 
Air Company, respectively, The catalyst used was 
Sherritt Gordon nickel powder C50-275, with 80°. of the 
particles having diameters between 100u and 150u. The 
superficial area of this powder is computed to be about 
55 cm? per gram. It had earlier been found “) that the 
catalytic activity of this powder was very low for some 
minutes at the first of each experiment. This effect, per- 
haps due to the presence of an oxide film, was found to 
be eliminated by pre-treatment in the ethylenediaminium 
carbonate solutions, under several atmospheres of hydro- 
gen pressure, and at the temperature of the subsequent 
expernnent, 

Experiments were carried out in a one-gallon auto 
clave, made by Autoclave Engineers Inc. Heating was 
done by a gas flame, controlled by a Wheelco controller, 
actuated by an iron-constantan thermocouple. Uhe tem 
perature was held constant to ©1°C. during each experi 


ment. Values were reproducible to ©1°C. Stirring was 
done by means of a stainless steel impeller operated con- 
tinuously at 700 r.pam, Uranium was added at the begin 


ning of each experiment, under operating pressure, in the 
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URANIUM CONCENTRATION 


0 10 20 30 40 50 60 
TIME (minutes) 


Figure 1—Typical reduction, showing two phases. 


form of a strong (30 gr. U/L) solution of uranium in 
ethylenediamine carbonate. Samples were withdrawn 
periodically, and analysed for dissolved uranium by means 
of x-ray fluorescence. These measurements were made on 
acidified solutions containing a thorium internal standard, 
using a Norelco x-ray fluorescence apparatus. The inten- 
sity of the uranium Ka line was compared with that of 
the Ka line from a known amount of added thorium. 
Calibration and reproducibility tests indicate that the 
analytical errors are 1% or less. 


Results and Discussion 


The course of a typical reduction is shown in Figure 
1, where two stages are clearly evident. The transition 
between the two stages occurred at various times and 
uranium concentrations, depending on the amount of 
catalyst present. The product was shown by x-ray diffrac- 
tion to be the same in both cases—a form of uranium 
diogide whose O:U ratio is somewhat greater than 2.0. 


Also evident from the data shown in Figure 1 is the 
independence of the second-stage reduction rate of uran- 
ium concentration. The linearity of the first stage is not 
evident in Figure 1 because of its being of such short 
duration. Rate measurements on the second phase reaction 
were made using the best straight line drawn visually 
through the experimental points. Although it is very difhi- 

cult to assess the accuracy of such measurements, the 
error in such a measurement is not likely to exceed 5% 
Reproducibility of experiments was not so good, however, 
and deviations of 10-15¢, were common. This can possibly 
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Figure 3— Dependence of reaction rate on hydrogen 
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Figure 2—Dependence of the second-stage rate on amount 
of catalyst. 


be attributed to variations in the activity of the catalyst 
surface, which is likely to be quite sensitive to the exact 
conditions of the experiment. 

The erent: of the reaction rate on catalyst area 
is shown in Figure 2. The evident first-order dependence 
is to be expected for a reaction which takes place entirely 
at the catalyst surface. 


The effect of hydrogen pressure on both the first and 
the second stages of the reaction is shown in Figure 3. 
The first-order dependence found for the first stage is 
in agreement with the results of Forward and Halpern. 
Two things are apparent regarding the second phase of 
the reaction: a marked increase in the scatter of the data, 
and a transition to a fractional order in hydrogen. Both 
these effects suggest that the reaction rate is controlled 
by some physical or mechanical process rather than by a 
chemical reaction. 


The rate of the reaction was studied as a function of 
temperature and an apparent activation of 8,700 calories 
per mole was indicated. Comparison of this value with 
that of 9,800 calories per mole found ‘*) for the first stage 
reaction suggests that the rate controlling processes may 
be similar. 

It was found that, in spite of the implication of 
Equation (1), the nature of the cation has a considerable 
effect on the rate of the reduction. As shown in Figure 4, 
the rates in ammonium carbonate and ethylenediaminium 
carbonate solutions are nearly the same, but about ten 
times faster than the rate in sodium carbonate solutions. 
The reason for this is not clear, but may lie in the fact 
that the faster rates are found in solutions which are 
strongly corrosive to nickel and maintain clean catalyst 
surfaces. 

In spite of the thick coating of uranium oxide, it was 
found possible to conduct several successive precipitations 
using the same batch of catalyst. In such cases, each new 
charge of uranium was added to the solution without 
interrupting the pressure, temperature, or agitation. 

Figure 5 shows a cross section of a catalyst particle 
coated with uranium oxide from five successive precipita- 
tions. Five layers are clearly visible, separated by darker 
regions. This darker material is senncianed with the con- 
tinued action of the hydrogen on the precipitate after the 
solution has been completely stripped of uranium. It is 
likely that this dark material is a more highly reduced 
form of the uranium dioxide. 
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Figure 4—Reduction rate curves for different cations. 


In Figure 6 is seen a cross section of a UO, layer 
deposited by two successive precipitations: the first of 3.0 
grams per liter, and the second of 6.0 grams per liter. 
From the observation that the thinner layer is on the 
outside, it may be concluded that the reduction proceeds 
at the original nickel surface, after transport of the re- 
actants through the deposited layer. It will be observed 
in Figure 6 that the outside layer is very thin compared 
to the inner one—much more than the expected factor 
of 2. This is not understood, although it may result partly 
from the fact that, being older, the outer lay er has been 
subjected for a longer time to the compressive effect of 
collisions between particles. 

Diffusion of reactants through this layer would sug- 
gest that the reaction rate should be dependent on the 
reciprocal of the total amount of uranium precipitated, 
in keeping with the well-known relation 

dn AC 


=D. be Peta Ieee = Ne , «€3) 
dt 6 ( 


where the thickness, §, of the diffusion layer depends on 
the mass of product precipitated. The reactant whose 
diffusion limits the rate is evidently molecular hydrogen, 
although the fractional order suggests that the control 
may be shared with some other reactant. 

If the density of the deposited material is constant, 
one would expect that the rate of the reaction could be 
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12 r 
Figure 5 — Photomicrograph of a polished section of a 


UO,-covered catalyst particle. (Cracking likely took place 
on drying). (X 720). 





Figure 6—Photomicrograph of a polished section of UO,- 
covered catalyst particles, showing the first-formed layer 


on the outside. (X 720). 


related to the amount, M, of uranium precipitated, by 
the equation 


d{U] — A 
dt ™M 


- 


where A is a constant. Figure 7 is a plot showing the 
rates of several successive reductions, plotted according 
to Equation (4). It is evident from Figure 7 that the data 
can best be fitted by an expression of the form 

a ee 5) 

=x = . 
The additional term, B, in Equation (5) is interpreted 
as representing the effect of abrasion or flaking of the 
UO, layer, so as to reduce the thickness of the diffusion 
barrier. Once a complete layer has been formed, this 
abrasion may be more or less constant. This postulated 
abrasive effect is supported experimentally by the obser- 
vation that an appreciable fraction of the uranium oxide 
is always found to be free of the nickel catalyst. 

The initial deposition of UO, appears to occur on 
active areas, such as grain boundaries and certain grain 
faces. This is shown in Figures 8 and 9, which are photo- 
micrographs of the surface of a nickel sheet used as 
catalyst for a very short time. 


(4) 
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uranium precipitated. 


From this it can be deduced that for a heterogeneous 
surface the transition from the first to the second stage 
of the reduction should be of the form 

log ee bt ae ...(6 
where [U], is the uranium concentration in solution at 
the point when the catalyst surface is completely covered 
with UO, This point can be identified experimentally as 
the first point where the data begin to fall on the straight 
line for the second stage, as indicated by an arrow in 
Figure 1. Experimental results support this deduction, as 


} 
can be seen from Figure 10, which shows the data from 


two separate experiments plotted according to Equation 
(6 

The gross features of the process seem now to be 
clear: The first deposition is at active sites, dependent 


on the surface structure of the nickel, giving a rate which 
is constant in tme. The deposition spreads laterally until 
the surface is covered, giving rise to a logarithmic transi- 
tion stage. Reaction proceeds by virtue of diffusion of the 
reactants through the deposited layer. This last step is 


cS 


altered in nature by a gradual flaking off or abrasion of 


I 


ch causes serious deviation of the 


the deposited layer, wh 


rate from the expected parabolic behavior. 





Figure 9—Vhotomicrograph showing the effect of grain 
boundaries on the catalytic activity of the surface. (X 16). 
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Figure 8—Photomicrograph showing the effect of crystal 
orientation on the catalytic activity of the surface. (X 16). 
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Nomenclature 

A 

B - = undefined constants. 

Cc 

AC = concentration difference across the diffusion barrier 

D = diffusion constant. 

AH = experimental activation energy. 

k = rate constant. 

M = mass of UO. precipitated. 

(dn/dt) = amount of material across the diffusion barrier per 
cm.” per sec. 

R. = rate of the second-stage reduction. 

6 = thickness of the diffusion layer 

[H»] hydrogen pressure. 

[Ni initial catalyst surface area. 

| uranium concentration in solution at time t. 

[I initial uranium concentration in solution. 

[ uranium concentration in solution at the point when 


the catalyst first becomes completely coated with 
uranium oxide. 
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Aeration in Tower-Type Fermenters' 


D. MURPHY?, D. S. CLARK® and C. P. LENTZ* 


The effects of air flow rate, oxygen partial pres- 
sure, aeration disk size and porosity, liquid height, 
and surface-active agents on aeration were studied 
in two unstirred towers or column-type vessels of 
different diameters using the sodium sulphite oxi- 
dation method. The coefficient of oxygen transfer 
increased with air flow rate up to the “flooding” 
point. The air flow rate required to cause “flooding” 
was increased by increasing the area of the aeration 
disk relative to the cross-sectional area of the tower. 
Within limits, aeration was improved by using disks 
of finer porosity. Under most conditions the coeffi- 
cient of oxygen transfer was much higher in the 
small tower. The coefficient of transfer was not 
affected by liquid height in the tower, but efficiency 
of absorption increased with height. Aeration was 
adversely affected by all surface active agents tested. 


ao was found to be the most important factor 
in scaling up a submerged fermentation of the beet 
sugar molasses for the production of citric acid in un- 
stirred tower- -type fermenters “!). The results of that study, 
in which the fermentation was transferred from a small 
glass tower (2.5 liters operating volume) to a larger one 
(36-40 liters operating volume), both aerated through 
sintered glass disks, indicated that more information on 
the physical factors governing aeration in this type of 
vessel should be obtained. 


Although a considerable amount of work has been 
published on mass transfer from bubbles leaving an orifice 
(2, 8. 4,5) on mass transfer in column-type gas-liquid con- 
tactors (67), on aeration in fermenters (*: 9% 1% 11,12), and 
in scaling up fermentation equipment (1% 1% 14, 15,16), jt is 
of indirect value in scaling up the citric acid fermentation 
mentioned above. Information directly applicable to the 
eqegenent and conditions in use should be obtained before 

1 further increase in scale is attempted. The work on 
mass transfer from bubbles, for example, shows that 
orifice size and shape affect bubble size and in turn mass 
transfer, but this information can not be used to predict 
the performance of a series of sintered glass disks. The 
work on tower-type gas-liquid contactors does not cover 
disk area in relation to column area, oxygen absorption, 
or a suitable range of disk porosities, and in addition 
most of the work is devoted to higher rates of gas flow 
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than are of interest in the fermentation under study. Most 
of the work on fermentation equipment deals with rela- 
tively short, mechanically agitated vessels rather than 
unstirred columns or towers (8 9 10, 11, 12, 13, 14, 15, 16), 

In the present work aeration studies were made with 
two sizes of glass towers similar to those used by Martin 
and Waters (7) and Steel et al). The work included 
determination of the effect of factors such as gas flow 
rate, disk area and porosity, column diameter and height, 
oxygen partial pressure, and antifoam agents on the 
absorption of oxygen by copper catalyzed aqueous sodium 
sulphite solutions in the towers (*). The sulphite oxida- 
tion method has been shown to provide a satisfactory 
indication of the relative effects of physical factors such 
as fermenter size, agitation, and air flow rate on the rate 
of oxygen absorption (12, 18, 18), 


Experimental 

The small tower consisted of a section of 2.5 in. (6.3 
cm.) diameter pyrex pipe 120 cm. in length with an addi- 
eal 6 cm. section containing the aeration disk at the 
bottom. The aeration section was made from a male 60/50 
pyrex joint fused to a Buchner funnel containing a 
sintered glass disk. A wet test gas meter was used to set 
the gas flow at the desired rate and a rotameter was used 
in maintaining a constant flow. The tower was filled to a 
depth of 96 cm. in most experiments, where greater 
depth was required, an addition section was added to the 
top of the tower. The temperature of the solution in the 
tower was thermostatically controlled at 25°C. during 
experiments by cool water flowing down the outside of 
the tower. Samples were drawn from the center of the 
tower through a glass tube set into the wall of the tower 
through a glass tube set into the wall of the tower about 
30 cm. above the aeration disk. 

The large tower was made from flange-connected 
sections of six inch (15.2 cm.) diameter pyrex pipe; the 
total height was 270 cm. The bottom section contained 

rubber covered steel plate in which cut-down Buchner 
funnels with sintered glass disks were set \!), This tower 
was filled to a depth of 192 cm. in most experiments 
Femperature was controlled at 25°C. by means of a U- 
shaped copper cooling coil in the tower. Sampling and 
gas flow regulating equipment was similar to that used 
with the small tower. 

Medium porosity disks of sintered pyrex were used 
in most experiments, With fine and coarse disks being used 
to determine the effect of porosity. Maximum pore sizes 
in the three grades of disk were 77.4, 12.2 and 5.4 microns 
9%), Except where the effect of disk area was being 
measured, one disk, 3.8 cm. in diameter and 4 mm. thick 
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Figure 1—Effect of gas flow rate on coefficient of oxygen 
transfer and efficiency of oxygen absorption. (Small tower 
with medium porosity disk covering 42% 
sectional area). 


of tower cross- 


] tower, and seven similar disks in the 
This represented about 45°, of the cross- 
A medium porosity disk 
large tower was also tested. 


was used in the smal] 
larger tower. 
sectional area of each 


the full diameter of the 


tower. 


The rate of oxygen transfer was determined by iodo- 
metric titration of the sulphite remaining in the reaction 
solution at various times during the test. Tap water was 
used to make up the reaction solution, and all chemicals 
were reagent grade. The initial concentration of sodium 
sulphite was 0.2 N in most tests, and the concentration of 
G U SO, 5H.,O was 8 x 10% M, well above the minimum 
level established by Fuller and Crist °°”). To prevent loss 
of sulphite, flasks and used in sampling were 


pipettes 


holed with nitrogen before use. Other measurements 
made included gas flow rate, pressure loss through the 


. a } 
aeration disk, and hold-up or 


volume of bubbles in the 
liquid column 
The cases velocity (S$.G.V., cm./sec.), the 


Jas 
J 
= 


coefficient of oxygen transfer (K,A, gm./(1.) (hr.) 
(atm.), and the “retention” (cm.*/1.) were calculated 
from the measurements made (see appendix). The 
superficial gas velocity was determined by dividing the 
volume rate of gas flow by the cross-sectional area of the 
towel! It was used instead of the volume rate of ¥as flow 
in order that results obtained with towers of different 
diameters could be compared directly. The coefficient of 


oxygen transfer was calculated by dividing 


Ox vee Nn absorbed pa 


the grams of 


hour per liter of sulphite solution 


UY thc LOY mean oF the iniet and outlet seateg pressures. 


This is < odified coefiicient which includes the surface 
area of the bubbles as well as the actual oak ient of 
transfer ‘*). The known partial pressure of the oxygen in 
tne yas enter g the tower and whe partial] pr ssure in the 
ge eaving the tower, cal ated from the sulphite ox) 
dation rate, were used in obtaining the Jog mean pressure 
(srap al integratior ed On pa wal pressure determina 
tions at vario ( eve the towers confirmed the 
alidity of y log mi pre ( alues at all except 
the highest and lowest superficial gas velocities. At thes« 
velocines the error was appreciable but not large enough 
to affect conclusions drawn from the results. The reten 


~ »® 
= 
< 
— 3 
od 
i 
20 
< 
- ° 
=x xe 
> 
© 
2 
wi 
oO 
uw 
uw 
Ww 





0.1 02 03 04 O58 O6 O7 


0.8 O09 1.0 


OXYGEN PARTIAL PRESSURE , ATM 


Figure 2—Effect of oxygen partial pressure on coefficient 

of oxygen transfer and efficiency of oxygen absorption. 

(Small tower with medium porosity disk covering 42% of 

tower cross-sectional area; superficial gas velocity 0.57 
cm./sec.). 


tion was determined by dividing the volume of gas held 
in the aerated column of solution by the number of liters 
of solution in the tower. The efficiency, or proportion of 
oxygen absorbed to that entering the tower, was also 
calculated for some tests. 


Results and discussion 
Superficial gas velocity and oxygen pressure 

In the small tower the coefficient of transfer (K,A) 
of oxygen from air increased linearly with the volume 
rate of gas flow until “flooding” occurred (Fi igure 1). In 
the flooding range there was a gradual change in bubble 
size and in the character of their movement; bubbles 
coalesced and became large and turbulence at the top of 
the liquid increased markedly. The predominance of large 
bubbles in the flooding range suggested that the decrease 
in K,A and hence efficiency was caused mainly by 
decreas bubble area. 


The coefficient of transfer decreased linearly with 
weep partial pressure of oxygen in the gas (Figure 

. Since approximately 80% of the oxygen entering the 
tower was absorbed for all partial pressure values tested, 
the decrease in K,,A can be ascribed mainly to a reduction 
in bubble area. However, the greater reduction in transfer 
area for bubbles of pure oxygen compared to bubbles with 
a lower oxygen partial pressure, for the same efficiency of 
oxygen absorption, explained these decreases in K, A for 
low — ial gas velocities only (Figure | The s rapid 
falling-off of K,A and efficiency at higher ish velocities 
with pure ae as compared to air suggested some 
other factor (Figure 1). 


Disk porosity and area 


The medium porosity disk produced much higher 
values of K,A (Figure 3) and retention (Figure 4) than 
disk. There was no appreciable difference, 
I<, A and retention values between the med- 
jum and fine disks. 
loss through the 


the coarse 
however, 
aeration disks 
sharply with decreasing pore size, 


Pressure increased 
but only slightly with 


rate of gas flow. At a superficial gas velocity of 1 cm./sec., 
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Figure 3—Effect of disk area and porosity on coefficient 
of oxygen transfer. (Small tower). 


for example, the pressure losses for the coarse, medium 
and fine disks were 8.5, 40 and 72 cm. of mercury, respec- 
tively, while the corresponding pressure losses at a 
velocity of 2 cm./sec. were 9, 44, and 76 cm. of mercury. 
The relativ ely small change in pressure loss with increas- 
ing rate of gas flow was caused by the increasing number 
of pores through which gas flow ed as the flow rate was 
increased. 

The main effect of increasing disk area relative to the 
cross-sectional area of the tower was to raise the flow 
rate at which flooding occurred (Figure The relation 
between disk size and the flow rate at which flooding 
started suggested that flooding was mainly a function of 
the gas flow per unit area of disk, although the rate at 
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Figure 5—Effeet of disk area and porosity on coefficient 
of oxygen transfer and efficiency. (Large tower). 
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Figure 4—Effect of disk area and porosity on gas reten- 
tion. (Small tower). 


which K,A decreased in the flooding ranges for the 
different sizes of disks indicated that the relative size of 
the tower was also involved. Retention increased approxi- 
mately linearly with superficial gas velocity for all disks 
over the entire range of velocities tested (Figure 4). 
Diameter and height of tower 

In the large tower (approximately twice the diameter 
of the smail tower) the effects of superficial gas velocity 
and of disk porosity and area on K,A, retention, and 
efficiency were similar to those in the small tower, but 
values of K,A under comparable conditions (except for 
tower diameter) below the flooding point were approxi- 
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Figure 6—Effect of disk area and porosity on gas renten- 
tion. (Large tower). 
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Figure 7 — Effect of sodium sulphite concentration on 
coefficient of oxygen transfer and gas rentention. (Small 
tower with medium porosity disk covering 42% of tower 
cross-sectional area; superficial gas velocity 0.24 cm./sec.). 


(compare Figures 3, 4, 5 and 6). In most instances, 
depending upon the superficial gas velocity and cross- 
sectional area of the disk used, the retention in the large 
tower was about 4 that in the small tower for comparable 
conditions. The superficial gas velocity at which flooding 
occurred was higher in the larger tower. With disks of 
small cross-sectional area (less than 40% of the tower 
cross-sectional area) the maximum values of K,A reached 
were about equal, but with disks the next size larger 
(Figure 3, 42; Figure 5, 45%) the maximum for the 
smal] tower was double the maximum for the large tower. 

Tests made with air over a range of liquid levels 
between 24 and 192 cm. in both large and small towers, 
showed that the values of K,A did not change appreci- 
ably with height if a correction were made for the effect 
of oxygen absorption on bubble area. Efficiency, of 
course, varied sharply with liquid height in these experi- 
ments because of the high rate of oxygen absorption by 
the sulphite solution. This effect is, therefore, a charac- 
teristic of the sulphite method. 


Sulphite concentration, antifoam agents, pulp suspensions 

Tests with initial concentration of sodium sulphite 
between 0.02N and 0.40N showed that the rate of oxida- 
tion was highly dependent on initial concentration in the 
range 0.1ON and 0.20N (Figure 7). In any given test the 
rate of oxidation remained constant until most of the 
sulphite was used, but this rate depended on the initial 
sulphite concentration. Tests in which sodium sulphate 
was added to the solution showed that the rate depended 
on the total amount of ~~. and sulphate rather than 
on sulphite alone (Figure , B, C), and the reaction 
rate was approximately the same for total concentrations 
over 0.20N. 

Photographs taken during tests using sodium sulphite 
solutions of different normality showed that bubble size, 
and hence total bubble area available for mass transfer, 
varied markedly with sulphite concentration (Figure 8). 
This probably accounted for most of the effect of sulphite 
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TABLE 1 
EFFECT OF SURFACE-ACTIVE AGENTS ON OXYGEN ABSORPTION 
(Tests made in small tower, medium porosity disk 42% of tower 
cross- s-sectional area in size, superfici ial | air velocity 0. 3 cm./sec.) 


| 





Agent | Type Concentration| KLA ‘enka 

(ppm) | (%) 

None _ - | 25 84 

Lard oil Antifoam 50 18 73 

| 400 18 73 

© Octadecanol | Antifoam 400 18 72 
in mineral oil 

Antifoam 60! Antifoam 10 6.0 34 

50 2.0 13 

100 eae 16 

400 2.4 15 

Aerosol OT? Wetting + 15 66 

agent 10 16 67 

Dreft$ Detergent 10 14 64 

20 12 58 











Senesel Electric Company, silicone “products Department, 
Waterford, N.Y., 

?American Cece Company, 
Division, New York, N.Y., U.S.A. 

3Procter and Gamble Limited, Toronto, Ont. 


Industrial Chemicals 


concentration on K,A. Bubbles did not appear to decrease 
further in size as the normality was increased above 0.175. 

All of the antifoam and surface-active agents tested 
affected aeration adversely (Table 1). The magnitude of 
the effect varied widely, depending on the agent used, 
but with the exception of antifoam 60, did not vary 
appreciably with the amount of the agent added. Under 
the conditions used these agents did not affect bubble 
size noticeably. 

The addition of small amounts of finely ground filter 
paper affected K,A markedly (Figure 9). The variation 
of K,A with gas flow rate with any given amount of 
pulp was similar to the relation under comparable con- 
ditions without pulp. 

Conclusions 

The results obtained in this study may be summarized 
as follows: 

(1) The coefficient of oxygen transfer (K,A)_ in- 
creased with gas flow rate up to the limit imposed by 
“flooding”. 

(2) K,A decreased with increasing oxygen pressure. 

(3) K,A increased with decreasing aeration disk pore 
size, within limits. 

(4) Aeration disk area relative to the cross-sectional 
area of the tower affected K,A markedly by raising the 
gas flow rate at which “flooding” occurred. 

(5) Values of K,A were generally much lower in the 
larger diameter tower than in the smaller one, 

(6) K,A was not affected appreciably by the height 
of the liquid column, but efficiency increased with height. 

(7) K,A was adversely affected by all of the surface- 
active agents tested. 

Comparison of the results obtained in the towers 
(maximum values of K,A 140-150) with those obtained 
in mechanically agit: mel vessels by Cooper, Fernstrom 
and Miller (12) ° maximum values of Ky equivalent to Ky A 

50-60) shows that the towers provide a relatively high 
degree of aeration. Nonetheless, the limits imposed by 
conclusions (1) to (6) above indicate that baffles, mechan- 
ical agitation, and elimination of antifoam agents may be 
necessary to obtain the required level of aeration for the 
submerged fermentation mentioned earlier “'), if it is to 
be successful in larger vessels. 
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APPENDIX 1 Q: oF 
‘ P Log mean partial pressure in liquid column 
2 Sample calculations : 
a ; ; P, P, 0.21 — 0.042 
The following calculations illustrate the methods used In (P3) - = 0.105 atmos. 
C determining SGV, K, A, efficiency and retention. ‘The particular P, 1.59 
example selected corresponds to points (0.87, 74.8) and (0.87, loge | 
- ) 


. 83.5) on Figure 1 and point (0.87, 135) on Figure 4. 


1. Known experimental conditions R 23.6 


‘Temperature of reaction SGV . oat a 0.87 cm./see. 
medium (1) 25°C. : * 
Partial pressure of Q, 7.85 ‘ 
oxygen 1n air entering K,A - - = 74.8 g 02/( 1.) (hr.) (atmos.) 
tower (P;) 0.2 | atmos. Ps 0. 105 
Volume of liquid in ne Q, 7.85 a al 
tower (V) +e. Iftheiency O. x 100 = rp x 100 = 83.5% 


Cross-sectional area of 





(EL) (A) (12.5) (27) 






tower (A) 27 sq. cm. : 
Airflow rate (R) 85.0 1./hr. (23.6 cm.*/sec.) Retention V >. 135.0 em.*/ 1, 
(Procedure as outlined in ‘“Procedure’’) e 
* * * 
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The Effect of Pressure on the Enthalpy of Gases’ 


REZA KORDBACHEN? and CHI TIEN* 


A generalized enthalpy-correction chart based on 
the law of corresponding states was constructed using 
the compressibility factor information reported by 
Nelson and Obert The range of pressures cov- 
ered the region from P, = 0.02 to P, = 40, and the 
temperature extended from T, = 0.265 to T, = 12.5. 
The results produced from the correlation of this 
work have been found to agree with those obtained 
from some of the more sophisticated methods already 
presented in the literature. 


HE pressure effect on enthalpy for any real gases has 
long been recognized. Watson and Nelson “) con- 
structed the enthalpy correction chart for hydrocarbon 
mixtures and subsequent work by Hougen and Watson ‘?) 
gave the first generalized chart for enthalpy correction. 


Their work is based on the principle of corresponding 
H* 

states and the numerical values of 1 

compressibility factor information reported by Dodge “). 

The inadequacy of this chart has been repeatedly indi- 

cated by the large deviation when compared with actual 

experimental dane. 


resulted from 


The recent contribution of Lydersen et al “) and also 
that of Curl and Pitzer ‘°) point to the inability of the law 
of corresponding states to treat all substances alike. To 
overcome this limitation they have suggested the intro- 
duction of a third parameter in addition to the conven- 
tional parameters of P, and T,. This is accomplished by 
using critical compressibility factor, Z,., in Lydersen’s 
work, and the so-called accentric factor in Curl and 
Pitzer’s work. Jn theory, this approach is undoubtedly 
superior to the conventional procedures already in exis- 
tence. However, the introduction of an additional para- 
meter increases the complexity of computations and 
requires additional information besides the critical tem- 
perature and pressure. In many cases, the third parameter 


is rather difficult to obtain. 
g calc ulations, there 


gener: alized chart based on two 


parameters. The purpose of this paper is to present such 


For certain preliminary engineerin 
still is a need for a simpler 


a chart capable of produci ing reliable enthalpy corrections 
PPP ETE EEE TEE E EEE EEE TET ETT TET Eee rssssss 


1Manuscript rece d December 20, 1958 

2Graduate student, Department of Chemical Engineering, Oklahoma State 
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sAssistant pro’ or, Department of Chemical Engineering, The University 
of Tulsa, Tulsa, Oklahoma, U.S.A 

Contribution from the Department of Chemical Engineering, The Uni- 
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without recourse to the use of the more refined approaches 
suggested recently. The greater accuracy is made possible 
by using the compressibility peat ‘information reported 
recently by Nelson and Obert 6 The superiority of 
Nelson and Obert’s work is aes by the fine agree- 
ment between the predicted values based on their work 
and the actual values. Also the work of fugacity coeffi- 
cient ‘*) has proved to be worthwhile. 


Development of Correlation 

The enthalpy, as a thermodynamic propery, can be 
considered as a function of temperature and pressure. 
The partial derivative of H (enthalpy) with respect to 
P (pressure) is given by 


OH ; _ fav 

; AE We i) is veces 

OP), OT)» 
Since at zero pressure all gases behave ideally, by 

integrating Equation (1), we have 


fav 
‘I ~ —V 2 ee ree (2) 
OT] p 1 


where H* and H are the enthalpy at temperature T for 
P =O and P respectively. 


p 
(H* — H)r = 


oO 


From the equation of state, we have 


ala Pee ed 
PV =ZRT or Z or ..(3) 
R' 


Differentiating equation (3) with respect to P and 
substituting into Equation (2), it gives 


rr 
nig? ( (9z\ ap, (4 
TT Nar), P, ) 


I.quation (4) served as the basis for the construction 
of the present chart. Numerical values of the compressi- 
bility factor z as reported by Nelson and Obert (7 were 


replotted as z vs T, with P,. as parameter. The partial 
OZ ‘ : 
) is then evaluated by graphical differ- 
Py ' 


H* H 


derivativ e,( 7 
OT, 


entiation using the technique described by Worthing and 
Geffner “), 
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Figure 


* 


. 07 j H { , 
With ( oT ) known, the quantity, 1 , can easily 
aT.) p. ’ 7 


be computed from Equation (4) by means of graphical 
integration. 


Results and Comparison 

Figures 1 and 2 are plots of the calculated results 
obtained in this work. These figures give the plots of the 
* of 

| 

sure, P,, with T, as parameters. The use of this general- 
ized chart is simple and straightforward. For any given 
compounds at specified temperature and pressure, it is 
only necessary to compute the reduced temperature and 


enthalpy correction term, . versus reduced pres- 


W* 
pressure and the subsequent value of —, can then be 


ci 
obtained from the chart. 


Comparisons were made between this work and some 
other existing methods which include those by Hougen 
and Watson ‘*), Lydersen et al ‘*) and by Curl and Pitzer 
‘5), Results of the preceding methods of computation were 
also compared with experimental data for carbon dioxide, 
methane, ethane, isobutane, n-pentane, and ammonia 

*, 10,11), These comparisons are presented in Table 1. 


164 


rl FAR 








AS 






AL wg 
YM 








Sy ve 
ern ven 





| LAWL 





| de 
lA 


7 8 9 10 15 20 25 30 40 
Reduced Pressure, P 
2. 


The comparisons indicate that in several cases the 
predicted value from this work is the lowest among all 
the existing methods. In the case of carbon dioxide. the 
agreement with experimental values is very good with 
an average deviation of less than 2%. This represents a 
better agreement than is possible from the work of Curl 
and Pitzer, Lydersen, Greenkorn and Hougen. The pre- 
vious two-parameter charts give an average deviation of 
approximately 7%. Similarly this work gives fair agree- 
ment in the case of ethane. The deviations for the cases 
of isobutane and n-pentane are somewhat greater, but it 
should be noted that the other methods also have cor- 
respondingly greater deviation for these two compounds. 


It would be of interest to compare results for two 
extreme cases of this method of computation with other 
methods: methane with Z, = 0.289, and ammonia with 
Z.. = 0.234. For methane, this work gives better agreement 
at low temperature while that by Hougen and Watson 
fits the experimental values better at higher temperatures. 
The values based on Lydersen, Greenkorn and Hougen 
give better agreement for T, = 2.068, but for lower 
temperatures this work gives better results. In all cases, 
the predicted values from this work as well as that by 
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TABLE 1 
COMPARISON OF PREDICTED AND EXPERIMENTAL VALUE OF (H* — H)/T, 


Compound: Carbon Dioxide P. = 73.0 atm. T. = 304.1°K Z. = 0.275 
Tr, = 1.1 Te = 42 
, : , , : 
% ( ( ¢ / c fe c ke ¢ i ( t 
P, | Expt. | C.P. | Dev.| H.W. | Dev. |L.G.H.) Dev. | K.1 Dev. | Expt.| C.P. | Dev. | H.W. | Dev. |L.G.H.| Dev. | K.T. 
0.5 | 0.953 1.05 }10.0 | 1.05 |10.0 | 6.93 |—2.0 | 0.98 3.0 | 0.775) 0.795} 2.5 | 0.88 |13.5 | 0.715|—7.6 | 0.772 
1.0 | 2.325 | 2.36 | 1.5 | 2.50 | 7.5 | 2.325, 0.0:| 2.36 1:5) 827 1 277 0.0} 1.90 | 7.5 | 1.73 |—2.2 | 1.73 
1.5 | 4.4 | 4.5 2.0 | 4.65 | 6.0 | 4.29 |—2.4 | 4.48 1.5 | 2.88 | 2.86 |—0.8 | 3.05 | 6.0 | 2.90 0.8 | 2.88 
2.0 | 6.14 | 6.24 | 1.5 | 6.50 | 6.0 | 5.72 |—6.7 | 6.20 1.0 | 4.11 | 4.09 |—0.5 | 4.20 | 2.0 | 4.01 |—2.5 | 4.05 
3.0 | 7.325 | 7.43 1.0 | 7.55 | 3.0 | 7.15 |—2.5 7.00 |—4.5 | 5.76 | 5.82 1.0 | 5.90 | 2.5 | 5.60 |—3.0 | 5.60 
| Av. +3.2] Av. +6.5| Av. —2.7 Av, +2.3 Av. +1.0| Av. +6.3} Av. —3.2/ Av 
| Ty. = 2.3 T, = 1.4 
cet es — 7 ee | - F | r = 7 : bine - i ; - an 
P, | Expt. cP. % % | % Ne % Te Te 
Dev. |H.W.| Dev. |L.G.H.} Dev. | K.T. | Dev. | Expt. | C.P. | Dev. | H.W. | Dev. |L.G.H.| Dev. | K.T. 
0.5 | 0.675] 0.655) —3 0 | 0.75)11.0 | 0.578) 14.2 | 0.65 |—3.5 | 0.576 | 0.536 7.10) 0.65 |13.0 (0.476 | 17.4 0.542 
1.0) £41 | 1.41 | @.0 | 1.550.060) 1.37 3.0 | 1.41 0.0 | 1.193 | 1.132 }—5.0 | 1.30 | 9.0 {1.112 |—7.0 |1.175 
1.3 | 2.25 | 2.19 |=2.0 | 2.4) 2.8 | 2.23 0.0 | 2.23 0.0 | 1.810 | 1.790 |-1.0 | 1.95 | 7.8 |1.790 |—1.0 |1.810 
2.0 | 3.06 | 3.1 Oe | O02 | 4.0 | oOo 0.0 | 3.00 | 0.05 |} 2.415 | 2.500 1.4 | 2.60 | 5.4 |2.385 |—3.2 |2.465 
$0. | 4.51 a2 | 29.0 | 4.7 4.5 | 4.47 |-1.0 | 4.51 0.0 | 3.595 | 3.595 0.0 | 3.90 | 8.4 |3.575 |—0.5 |3.700 
Av. +6.8 | Av. +7.6 Av. +3.6 Av. —0.8 Av. 3.0 | Av. +8.8 Av. +5.8 Av. 
Compound: Methane P. = 45.8 atm. Fe. = 443.3°R Z. =,0.29 
Tr = 1.54 T, = 1.893 
P, Expt. |L.G.H. % Dev. | H.W. | % Dev.| K.T. | % Dev. | Expt. | L.G.H. | % Dev.| H.W. | % Dev.| K.T. | % 
0.297 | 0.238 | 0.28 +17.6 | 0.32 +34.5 0.22 7.6 0.163 0.117 —28.2 0.22 +24.5 
0.891 0.853 | 0.70 —17.9| 0.96 +12.5 0.83 nua 0.541 0.415 —23.3 0.64 +18.3 0.38 — 
1.485 1.589 | 1.40 —11.9 1.57 + 1.2 1.42 —10.6 | 0.904 0.750 —17.0 | 0.99 + 9.5} 0.72) — 
2.240 | 2.208 | 2.10 — 8.0; 2.40 + 5§.1 2.20 — 3.6 1.328 1.150 —13.4 1.45 + 9.6) 1.13} -— 
4.458 | 3.858 | 3.45 —10.6; 4.00 + 3.7 3.80 - 1.5 2.395 2.010 —14.0 |} 2.40 + 0.2} 2.2 
Av. +13.2 Av. +11.4 Av. + 5.2 Av. —16.1 Av. +12.4 Av. - 
T, = 2.068 
| C C Oo 
. \ - 7. 
P, | Expt. |L.G.H.| Dev. | H.W.| Dev. K.T. | Dev. 
0.297| 0.144 | 0.10 | —30.6 | 0.19 | +31.9 
0.891 | 0.452 | 0.38 15.9 | 0.55 +21.7 | 0.25 | 44.7 
1.485 | 0.756 | 0.62 | —16.5 | 0.86 | +15.3 | 0.52 | 30.3 
2.228} 1.095 | 0.92 —16.0 | 1.25 +14.2 | 0.90 | 17.8 
4 all 1.976 | 1.80 | — 8.9 | 2.00 | + 1.2 | 1.70 | 14.0 
| Av. —17.6 Av. +17.0 | Av. +27.0 
Compound: N-Pentane P. = 33.3 atm. rT. = 469.8°K Z, = 0.269 
| lr, = 1.088 
| 
Py Expt. | C.P. % Dev. H.W. o Dev. L.G.H. % Dev. Kk. % D 
2.04 6.95 6.85 1.4 6.82 1.9 6.21 10.6 6.67 2 
3.06 | 7.85 7.44 ~§.3 7.78 0.9 7.37 6.1 7.17 8 
4.08 8.19 7.94 3.1 8.02 ait tute 5.4 7.39 9 
6.12 8.47 8.32 —1.8 8.32 1.8 8.05 5.0 7.56 10 
10.20 8.52 8.5 0.3 8.01 6.0 7.45 12 
16.32 8.11 7.41 8.6 7.42 8.5 6.83 15 
| | 
| 
| Av. —2.6 \v. —2.6 Av. 7.9 Av. —10 
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Table 1 (Continued) 


Compound: Ethane P, = 48. ‘ T. = 305.2°K 


Z, = 0.285 
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Compound: Isobutane 36.0 atm. = 408.1°K 
| 
| 


OQ 


Expt. 


0.795 
.590 
. 880 
830 
.440 
335 
665 
850 


| 
tN 


| 
Owrneee OO | 
SNS OO FS ROS © 
SSOWWWEDN | 
NVMON WNW | = 


NSA We 
COCOONTA Uwe © 


p 


Compound: Ammonia 


P, = 111.3 atm. 


P, | Expt.* |L.G.H. 


0.0611 | 0.198 | 0.23 b | 3-0 | 0.155 : 
0.0917 | 0.295 | 0.326 i : we 0.250 | —15.: 
0.1222 | 0.398 | 0.46 a : -5 | @0.3552 | =i. 
0.1834 | 0.612 | 0.73 | oo .625 0.560 | — 8.: 


| | Av. +15. Av. +3.8| Av. —14.3 


*Values of H* are obtained from linear extrapolation with the information in 
Tables of Thermodynamic Properties of Ammonia, U.S. Bureau of Standards, 
Circ. 142, 1923. 


Lydersen, Greenkorn and Hougen are 10-15% lower than 
the experimental value, while those by Hougen and 
Watson are 10-15, greater. For the case of ammonia, the 
two-parameter chart of Hougen and Watson appears to 
be superior. However, it should be mentioned that the 
enthalpy of ammonia at ideal state (P, = 0) is obtained 
by linear extrapolation. Conceivably, this will cause certain 


; ; { ; . : 
error in the expression , Which might distort the 


T 
true picture. 

A question which arises in the case of methane is why 
a more refined method such as the one used by L ydersen, 
Greenkorn and Hougen does not give any better results 
than the simpler two-parameter method. It might be said 
that the experimental values used for comparison are not 
too accurately determined. If this is the case, some new 
experimental investigation should be in order. On the 
other hand, if the experimental values are reliable, the 
only conclusion would be that Z, might not be an ideal 
parameter to characterize the intermolecular effect, a 
comment which has been voiced by some current critics. 


166 


Oke BONN 

e e t 

Auk PwNn— 
Sw 
HAurwn 


Sweat 
SS 
OME Ee ONN ee 


an 
a 
w 
= 





LG.H. 


= 
— 
Y 
@ 
< 


| 
| 
| 
| 
| 
| 


1 | 
rm | 


| 
NON OUA~IWh 


62 
.53 
.80 
.59 
.90 
.15 
.65 
.02 


| 


eNO COCO 
WOOUNMNOUwW 


a 
to 


MUAMNI OWS 
on Ne © 

Pm OMNwWOoO NH 7 
NSNINTOA We © | 


NR 
© 
° 


Conclusion 


From the comparison of this work, the present inves- 
tigation makes possible the better prediction of enthalpy 
of gas from ideal behavior than other available methods 
and particularly at low temperatures. The work by 
Hougen and Watson appears to be superior only in the 
case of ammonia. 


Nomenclature 


H 
H* 
P 
P, 
P, 


Enthalpy of real gas 
Enthalpy of ideal gas 
Pressure 

Critical pressure 
Reduced pressure 

Gas law constant 
Temperature 

Critical temperature 
Reduced temperature 
Volume 
Compressibility factor 
= Critical compressibility factor 


uu dud 


ou uw il 


References 


(1) Wason, L. M., and Nelson, E. F., Ind. Eng. Chem., 25, 880 

(1930). 
Hougen, O. A., and Watson, K. M., Chemical Process Principles, 
Part 2, John Wiley and Sons, Inc., New York. 

) Dodge, B. F., Ind. Eng. Chem., 24, 1353 (1952). 
Lydersen, A. L., Greenkorn, R. A., and Hougen, O. A., Generalized 
Properties of Pure Fluids, Report No. 4, Engineering Experiment 
—. a of Wisconsin, October, 1955. 
Curl, and Pitzer, K. S., Ind. Eng. Chem., 50, 365 (1958). 

3) ea on C., and Obert, E. F., Tans. ASME, 76, 1057 (1954). 
Nelson, L. C., and Obert, E. F.. AIChE. Journal, 1, 75 (1955). 
Honari, S., and Brown, G. G., Private Communication. 
Worthing, A. G., and Geffner, J, Treatment of Experimental Data, 
John Wiley and Sons, Inc., New York, 1943. 
ao, B. H., Olds, R. H., and Lacey, W. N., California Oil World, 
(1946). 

) Tables of Thermodynamic Properties of Ammonia, U.S. Bur. Stand- 
ards Cire, 142, 1923. 

*x* && ® 


The Canadian Journal of Chemical Engineering, August, 1959 








